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A STUDY OF ROTOR BROADBAND NOISE MECHANISMS 
AND HELICOPTER TAIL ROTOR NOISE 

Shau-Tak Rudy Chou, Ph . D . 

Cornell University, 1987 


A study is made of the relative problems of rotor broadband noise 
mechanisms and the helicopter tail rotor noise. The rotor broadband 
noise mechanisms considered in the present study are (1) lift 
fluctuation due to turbulence ingestion, (2) boundary layer/trailing 
edge interaction, (3) tip vortex formation, and (4) turbulent vortex 
shedding from blunt trailing edge. Predictions are compared to 
available experimental data and show good agreement. The study shows 
that inflow turbulence is the most important broadband noise source 
for typical helicopters' main rotors at low- and mid- frequencies . 
Trailing edge noise and tip vortex noise are found to be important at 
high frequencies; they are very sensitive to the rotor blades' angle 
of attack. Trailing edge thickness noise is also very important; it 
generates a large spectrum hump, and is very sensitive to any change 

of the trailing edge thickness. 

Due to the size difference, isolated helicopter tail rotor 
broadband noise is not important compared to the much louder main 
rotor broadband noise. However, the inflow turbulence noise from a 
tail rotor can be very significant because it is operating in a highly 
turbulent environment, ingesting wakes from upstream components of the 


xi 


helicopter. Our study indicates that the main rotor turbulent wake is 
the most important source of tail rotor broadband noise. Tail rotor 
harmonic noise due to nonuniform inflow resulting from various 
upstream disturbances (e.g. main rotor, fuselage, engine exhaust, and 
main rotor hub) is also investigated. The harmonic noise due to the 
main rotor mean wake is found to be very important. 

Finally, the tail rotor harmonic noise due to ingestion of main 
rotor tip vortices is studied. The CAMRAD code is used to find the 
main rotor tip vortex trajectories; then, using the near-normal 
blade-vortex interaction model of Amiet, the noise radiation is 
calculated. Tail rotor position is found to be very important as it 
directly affects the geometry of the interactions. The tail rotor 

phasing is also found to be important in case of commensurate tail and 
main rotor RPM's. 
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Chapter I 


INTRODUCTION 


In those days, the world teemed; the people multiplied, 
the world bellowed like a wild bull, and the great god was 
aroused by the clamour. Enlil heard the clamour and he said 
to the gods in council, "The uproar of mankind is 
intolerable and sleep is no longer possible by reason of the 
babble." So the gods agreed to exterminate mankind. (From 
a translation by N. K. Sandars of The Epic of Gilgamesh , a 
Babylonian story believed to be the forerunner of the flood 
story in the Bible.) 


In recent years, helicopters have been proved to be a convenient 
and reliable means of transportation. Their ability to land, take 
off, and maneuver in areas inaccessible to fixed wing airplanes has 
been successfully used in rescue work and inter-city commuter 
applications. However, the increasing use of helicopters in both 
civilian and military applications has raised attention to the noise 
they generate. In civilian applications, the high noise level has 
prejudiced the very operation which helicopters are uniquely fitted to 
perform - inter-city transportation. In military applications, the 
far- field noise provides unnecessary early warning of the approach of 
helicopters to their potential enemies. The internal noise problems 
are also significant, causing helicopters to be unattractive to 
prospective civilian passengers and, in the military field, internal 
noise levels often substantially exceed the acceptable limits (Lowson 
and Ollerhead, 1969). 
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In this thesis, chapter II will review and derive basic governing 
equations of aeroacoustics . Chapter III will investigate mechanisms 
that generate rotor broadband noise; various analyses will be extended 
or developed and they will be shown to be capable of predicting rotor 
noise spectra. Chapter IV will study the harmonic and broadband noise 
from helicopter tail rotors using existing aerodynamic and acoustic 
analyses. Chapter V will examine the interaction between helicopter 
tail rotor and main rotor tip vortices and the noise it generates. 

The primary noise sources for a typical helicopter are the main 
and tail rotors, engine(s) , and gearbox. According to the origins of 
these noise sources, they can be divided into those generated 
aerodynamically (i.e. noise from main and tail rotors) and those 
generated mechanically (engine and gearbox noise) . Cox and Lynn 
(1962) found that at moderate distance from a typical helicopter, the 
various noise sources can be listed according to their order of 
importance : 

(1) blade slap (when it occurs) 

(2) tail rotor rotational noise 

(3) main rotor broadband noise 

(4) main rotor rotational noise 

(5) gearbox noise 

(6) turbine engine noise 

(7) other sources 

In general, the relative importance of these mechanisms depends 
upon the particular helicopter design and the operating conditions. 



The main rotor impulsive noise (i.e. blade slap) is the most important 
helicopter noise source when it occurs. However, in the case of 
modern helicopters, main rotor impulsive noise sometimes is not as 
significant due to careful design. This often leaves the rotational 
noise from the tail rotor and the broadband noise from main and tail 
rotors as the most important noise sources on helicopters. In the 
present study, we will focus our attention on these particular 
aerodynamic noise sources. 



Frequency (Hz) 

Figure 1 . 1 : Filter Characteristics for A-Weighted Sound Level. 



-4 - 


The frequencies of interest in the present study are mainly 
determined by human annoyance (or detection in some cases). The 
common measures of annoyance, such as the perceived noise level (PNdB) 
or A-weighted sound level (dBA) , account for the fact that humans find 
low frequencies, say below a few hundred Hertz, much less annoying. 

The filter characteristics for A-weighted sound level are shown in 
Figure 1.1. On the other hand, if long distance propagation is a 
factor for the rotor in question, then high frequencies can be 
attenuated significantly by molecular absorption. For example, after 
propagation over one kilometer, sound waves with frequencies above a 
few thousand Hertz are attenuated drastically (George, 1978). Thus 
frequencies in the range from approximately one hundred to several 
thousand Hertz are of primary interest. 

Figure 1.2 shows a typical acoustic time history of a helicopter 
in flight. Two distinctive signal patterns are apparent: the first 

is periodic impulses with fundamental frequency originating from the 
blade passing frequency, and the second is essentially a random 
background signal. After Fourier analysis, the difference between two 
types of signals becomes even wider. The periodic impulses generate a 
line spectrum at the fundamental frequency (i.e. blade passing 
frequency) and its harmonics; the random signal leads to a continuous 
but possibly peaked broadband spectrum, A typical helicopter noise 
spectrum is given in Figure 1.3. Generally speaking, the periodic 
impulses are more distinct during high speed forward flight and are 
due to the effects of the blade steady loading, due to volume 
displacement effects, due to the periodic blade-vortex interactions, 
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FULL SCALE 




, V 30° 



FULL SCALE-UNAVERAGED 



Figure 1.2: Measured Pressure -Time History for a UH-1H in 60 Knot 

Partial - Power Descent (400 ft/min Rate of 
Descent), <f> - -30°, from Schmitz and Yu (1983). 


and due to other periodic blade loading variations. On the other 
hand, the random background signal is generally due to various 
blade - turbulence interaction effects. They can be the primary noise 
sources in the absence of extensive impulsive noise, as in the hover 
case. A series of hover spectra are given in Figure 1.4. 

Despite extensive research over the past fifty years, and 
particularly over the last fifteen years, the relative importances of 
various rotor aerodynamic noise mechanisms are only now becoming 
understood. The accuracies of the existing analyses are also hard to 
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INDICATED AIRSPEED = 80 knots 
RATE OF DESCENT = 0 ft/min 



Figure 1.3: In- Plane UH-1H Acoustic Power Spectrum for 80 Knot Level 
Flight, from Schmitz and Yu (1983). 


document. The primary reason for these difficulties is that there are 
a large number of noise mechanisms on rotors which can be important in 
different parts of the acoustic spectrum, depending upon the rotors' 
parameters and their operating environment. The wide variety of 
source mechanisms is due to various aeroacoustic effects: boundary 

layers, separated flow, and inflow turbulence; high Mach numbers, 
including nonlinear effects; blade-vortex interactions; non-uniform 
inflow, etc. (George, 1978). These mechanisms have been postulated 
and studied both analytically and experimentally to access their 
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Low frequency broodborxj noise High frequency broodband noise 

(20 Hz bandwidth) ( 100 Hz bandwidth) 

Lift 3050 lb, 160 rotor rev/mm 



Figure 1.4: Broadband Noise Spectra, Full Scale S-55 Rotor, <j> * -75 , 

from Leverton (1973). 

possible significance as noise sources. In general, the mechanisms 
each affect different parts of the acoustic spectrum. Moreover, on 
craft with either tandem or main and tail rotors, many of these 
mechanisms can interact with each other and between rotors. Thus, in 
many cases, it is not clear which mechanisms are dominant in many 
operating conditions for full scale helicopters, propellers, etc. A 
brief summary of the various noise generating mechanisms will be given 
in the remainder of this chapter where emphasis will be put on 
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explainmg the physical phenomena rather than on deriving the theory. 
The mechanisms will be analyzed in the later chapters. 

Discrete frequency noise (sometimes called rotational or harmonic 
noise) is caused by steady or harmonically varying forces, volume 
displacements, or nonuniform inflows. At low to moderate blade tip 
Mach number these can be due to the basic blade rotation and forward 
flight of a helicopter or to the steady inflow variations. These 
mechanisms have been analyzed by a number of researchers (Gutin, 1936; 
Deming, 1938; Hubbard, 1953; Lowson, 1969; and Wright, 1969). 

The first theoretical model of rotor noise was developed by Gutin 
(1937) who recognized that steady aerodynamic forces on a propeller 
blade act as acoustic dipoles. His analysis showed that the first few 
noise harmonics are related to the steady thrust and drag (torque) 
forces acting on a propeller blade. 

For many years Gutin noise" formed the sole basis for rotor 
noise prediction. Since the discrete noise spectrum due to steady 
loadings decays very rapidly with frequency, the steady loading noise 
is generally restricted to the first dozen or so harmonics of the 
blade passing frequency. Thus it is not usually of importance to 
helicopter main rotors or large wind turbines as these frequencies 
generally lie in the frequency range below 100 Hertz where human ears 
are not very sensitive. These low order harmonics are, however, very 
important to high speed propellers or tail rotors. 

Since the rapidly decaying steady loading noise severely 
underestimates the harmonics at high frequencies, it was felt that 
sizable unsteady loading fluctuations due to steady but spatially 
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nonuniform inflow to the rotor generate noise at high frequencies. 

This problem was formulated and solved by Lowson and Ollerhead (1969) 
and Wright (1969) who analyzed the noise radiation due to azimuthally 
varying blade loadings which are steady in time. They found that the 
higher loading harmonics of rotor blades are extremely important to 
the high frequency discrete frequency noise. In fact, at high 
frequencies the sound generated by even very small loading harmonics 
dominates that generated by the steady loading noise analyzed by 
Gutin. 

Impulsive noise (sometimes called blade slap) consists of nearly 
distinct repeated pulses at blade passing frequency. It is a special 
type of rotational noise. When an impulsive acoustic time history is 
Fourier analyzed, the repeated pulses will yield discrete or harmonic 
spectra, but their particular identity is due to their impulsive time 
histories and origins. These pulses are caused by particular events 
at certain blade azimuth angles such as main rotor blade -vortex 
interactions or high speed noise due to local transonic blade motion 
toward the observer (say Mach number greater than approximately 0.75). 
These noise mechanisms have been analyzed by Widnall (1971), Ffowcs 
Williams and Hawkings (1969a), Farassat (1975, 1980, 1983, 1984), 
Hanson (1979), Yu and Schmitz (1980), George and Chang (1983c, 1984b), 
and George and Lyrintzis (1986). Impulsive noise is unquestionably 
the most important noise source on helicopters or wind turbines when 
it exists. However, a prime goal of aeroacoustic rotor design or 
operation is to avoid impulsive noise generation by controlling the 
blade -vortex interactions and by avoiding high tip 


Mach number. This 


- 10 - 


often leaves broadband noise as the important controlling noise in 
many situations where relative tip speeds are low and blade-vortex 
interactions are not intense. 

Broadband noise is characterized by a continuous (although 
sometimes humped or peaked) spectrum and is caused by disturbances 
which are not precisely repeated at each blade revolution but are 
basically random in nature. These random disturbances are generally 
due to some sort of turbulence/rotor blade interactions. The 
turbulence can either be that already existing in the surrounding 
atmosphere or can be generated by the blades' motions. Recent reviews 
of broadband noise research may be found in the works of George and 
Chou (1983a, 1984a) and Schlinker and Brooks (1982). 

Opinions regarding the origin of broadband noise vary somewhat 
among investigators. Part of the reason for this variation is that 
broadband noise actually has several origins, although they are all 
related to the random loading/surface pressure fluctuations on rotor 
blade due to interactions with turbulence. Several possible causes 
have been proposed and investigated. Figures 1.5 and 1.6 give the 
sketches of the aerodynamic aspects of some of the mechanisms. 

One of the sources of the random part of the rotor broadband 
noise is the fluctuating loading associated with the encounter of the 
ambient turbulence (Sharland, 1964). When a rotor blade is moving in 
turbulent flow, the local incidence angle is determined by the 
resultant of the mean velocity and the instantaneous turbulent 
velocity. These random fluctuations in incidence angle lead to 
unsteady loading fluctuations on the rotor blade and subsequently 
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Figure 1.5: Broadband Noise Sources of a Stationary Airfoil . 

generate noise. This mechanism is particularly important in the low- 
to mid- frequency range. The incident turbulence may be that already 
existing in the atmosphere, or it may be generated by the 
recirculation wake due to ground effect, or it may be from the 
turbulent wake of the preceding blades. 

Another type of random force fluctuations on surfaces in a moving 
fluid is the surface pressure field arising from an attached turbulent 
boundary layer (Sharland, 1964). The problem regarding the sound 
generation from boundary layer turbulence was solved by Powell (1959). 
He used the "reflection principle" to show that the major surface 
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dipoles vanish on an infinite, flat, and rigid surface. This leaves 
only the viscous dipoles with their axes lying on the surface remain 
effective. Since such viscous stresses can only become significant at 
low Reynolds numbers, direct radiation from the turbulent boundary 
layer is generally not important. However, Hayden (1972) suggested 
that, providing the acoustic wavelength is much smaller than the blade 
chord, pressure fluctuations in the turbulent boundary layer can 
interact with a sharp trailing edge to produce significant noise 


radiation . 
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Force fluctuations can also appear when there is vorticity shed 
from the rear of a moving body. A classical example of this is the 
"Karman "Vortex Street" occuring behind a circular cylinder in certain 
Reynolds number ranges. Although the shape of rotor blades is 
different from that of a circular cylinder, similar pressure 
fluctuations which are associated with nearly periodic vortex shedding 
occuring in the laminar boundary layer Reynolds number range. The 
nearly periodic nature of the force fluctuations gives rise to a 
continuous but peaked spectrum, this is often identified as "vortex 
noise" or "high-frequency broadband noise" (Paterson, et al . , 1973). 
However, this phenomenon generally does not exist in full-scale rotors 
as the flows around the blades are usually turbulent. 

The importance of turbulence in blade tip flows was also 
suggested by several investigators. Lowson's experiments (1972) on 
low speed axial flow fans indicated that tip shape has a significant 
influence on the high frequency broadband noise. He suggested that 
this can be traced to the breakdown of the tip vortex, which induces 
high intensity turbulent flow over the outer portion of the rotor 
blades . 

A similar mechanism of turbulent vortex shedding from blunt 
trailing edges has been identified by Brooks and Hodgson (1980) . 

Unlike lamiflar vortex shedding, the vortices shed from the trailing 
edges are essentially random, although some correlation exists between 
the pressure fluctuations and the free stream velocity and the 
thickness of the trailing edges. 
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Other sources of broadband noise from rotor blades can be due to 
turbulence in locally stalled regions (Paterson, et al . , 1975). This 
phenomenon can be associated with either high angle of attack or close 
encounter with a vortex shed from previous blade. Presently, there is 
no analytical model available for either the local separation or the 
noise radiation. 

Most of the discussions above are concerned with the noise 
generated from helicopter main rotors. Despite its smaller size 
compared to the main rotor, the tail rotor's importance on the overall 
noise should not be underestimated. It has been known for some time 
that the tail rotor is an important source of helicopter noise. 
Although perhaps an exaggeration, Lynn et al. (1970) reported: "For 
nearly all flight conditions the tail rotor is the predominant noise 
source for single rotor helicopters." Similarly, Leverton (1977, 
1980), Balcerak (1976), Levine (1976), and Pegg and Shidler (1978), in 
their analytical and experimental studies on tail rotor noise, have 
all emphasized the importance of tail rotor noise due to interaction 
with the wakes of main rotor and other upstream objects. 

Tail rotor rotational noise is of particular importance since 
tail rotors are often operated at higher speeds than main rotors . The 
rotational noise of a typical helicopter tail rotor lies in the 100 to 
1,000 Hertz range, a range which is extremely important to audibility 
and annoyance. Tail rotor broadband noise is also very important due 
to the highly turbulent inflow which originates from various 
aerodynamic wake effects. 
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As mentioned previously, the added difficulty involved in the 
study of tail rotor noise is the complexity of the inflow. The inflow 
field of the tail rotor is very difficult to predict, it consists of a 
number of components: the overall main rotor wake, the main rotor tip 

vortices, the wakes of fuselage and main rotor hub, and the engine 
exhaust flow. These components affect both the mean flow and the 
turbulence properties as seen by the tail rotors. The overall main 
rotor mean wake can be calculated by several methods but it has a 
fairly minor influence on the tail rotor noise, because it leads 
primarily to low order loading noise harmonics. On the other hand, 
the concentrated main rotor tip vortices are quite localized and will 
lead to more annoying, higher frequency noise. The trajectories of 
the main rotor tip vortices are difficult to predict; they follow 
initially helicoidal paths which are perturbed by: (1) the roll-up of 

the vortex wake into a horseshoe vortex system; (2) the self- induced 
instabilities of the vortex trajectories; (3) and the action of the 
pre-existing atmospheric turbulence. The actual structures of the 
vortices are also not known except that they are generally turbulent. 
They can contain axial velocities in their cores which are very 
important to the unsteady loading fluctuations on the tail rotor 
blades. Also under certain conditions the main rotor tip vortices are 
found to be' "diffuse" or "burst" (the "vortex breakdown" phenomenon). 
When this occurs, the vortices become more spread out and turbulent. 
This strongly affects the characteristics of both the mean and 
turbulent inflow seen by the tail rotor. 
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In the present study, we limit our attention to the aerodynamic 
noise generation from helicopter rotors. In particular, the broadband 
noise from main and tail rotors, and the harmonic noise from tail 

rotors due to wakes from main rotor, hub, fuselage, and engine exhaust 
are studied thoroughly. 

Chapter II treats the general problem of aerodynamic noise 
generation. Various forms of governing equations are derived and 
their solutions discussed. For the simplified cases, the sound field 
from moving acoustic point singularities (in both straight-line and 
circular motion) is discussed. Finally, several approaches for 
predicting rotor noise are reviewed to conclude this chapter. 

Chapter III addresses part of the helicopter noise problem - 
broadband noise, particularly those from helicopter main rotors. It 
reviews previous broadband noise analyses and describes the extensions 
we make in the present study. For inflow turbulence noise, we extend 
the analysis of George and Kim (1976) to allow the use of von Karman 
spectrum. For boundary layer/trailing edge noise, the present study 
extend the analysis of Kim and George (1980) to include the effect of 
rotor blades' angle of attack. We also extend the analysis of George 
et al. (1980) for the tip vortex formation noise; a 3-dimensional 
normalized surface pressure spectrum is obtained. We also study the 
noise due to turbulent vortex- shedding from blunt trailing edges and 
develop an analysis to predict it analytically. Calculations based on 
various analyses of broadband noise mechanisms are compared to each 
other and to available experiments. The aims of this chapter are to 
help understand which broadband noise mechanisms are important under 
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which circumstances, to identify a number of satisfactory, existing, 
and well - documented experimental measurements, and to evaluate the 
various analytical approaches. It is shown that several satisfactory 
analytical approaches are available and their limitations are 
delineated. These analyses can show which mechanisms are important in 
which cases and are able to predict absolute spectra to within about 5 
dB for clean experiments. 

Chapter IV describes results from a comparative study of 
helicopter tail rotor noise. With components including the wakes of 
the main rotor, the main rotor hub, the main rotor shaft, the fuselage 
boundary layer, the fuselage separation, and the engine exhaust, tail 
rotor inflow is very complex and turbulent. This leads to higher than 
normal noise radiation considering the tail rotor's size. The primary 
harmonic noise source we focus on in this chapter is the interaction 
between the tail rotor and the main rotor mean wake. For the 
broadband noise, the intensities and length scales for turbulence 
generated by various upstream disturbances are estimated and the 
resulting noise from a typical tail rotor is calculated and discussed. 

In chapter V, we focus on the harmonic noise generated from tail 
rotor blades chopping the main rotor tip vortices. The analysis of 
Amiet (1984, 1986a) for near-normal blade-vortex interaction is 
reviewed. The details of a complicated simulation we devised for the 
present study are outlined. Results of the noise calculations for 
full-scale and model-scale UH-1 helicopters are discussed to conclude 

this chapter. 

The summary and conclusions drawn from the present research are 


given in Chapter VI . 


Chapter II 

BASIC AEROACOUSTIC FORMULATION 


This chapter reviews the basic theory of aerodynamic noise 
generation. Emphasis is placed on the mathematic modelling of the 
sound fields generated by given noise sources rather than on the 
physical modelling of the aerodynamic noise sources, which will be 
discussed in later chapters. Various forms of the governing equations 
will be derived and the general solutions will also be discussed. To 
conclude this chapter, several approaches developed previously to 
appiy the basic acoustic equations to rotor noise prediction are 
reviewed and discussed. 

2.1 CLASSICAL THEORY 

The theory of sound received its fair share of attention during 
the advance of modern science. The fundamentals of acoustic theory 
had been set down well before the end of nineteenth century (Rayleigh, 
1896). These fundamentals were applied to the studies of a number of 
harmonic phenomena, e.g. tuning forks, organ pipes, and church bells, 
etc. Basically, the main interests of the classical acoustic theory 
were the production and propagation of small disturbances in the 
pressure (or density) caused by stationary sources in an unbounded, 
stationary medium. Rayleigh's basic theory has been used by a number 
of researchers through a wide span of time. This not only gives 
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witness to the correctness of the works of Rayleigh but also indicates 
that acoustics has, in its fundamentals, changed far less since 
Rayleigh's time than other disciplines in science. 

For the basic theory of aerodynamic sound, the main facts that 
needed to be recalled from the classical acoustic theory are the 
elementary solutions of the equation of motion, namely, sources 
(monopoles), and dipoles (Rayleigh, 1896; Morse, 1948). The classical 
theory of acoustics correctly identified that sound can be generated 
by the injections of either mass or momentum into the acoustic medium. 

For a point monopole (injection of mass), it is the variation of 
mass outflow from the source that generates sound (in the present 
study, we generally call source of this type monopoles rather than 
sources as in classical acoustics). Assume new fluid is introduced at 
a rate Q(t), we may call its time derivative Q(t) the strength of the 
monopole, since the pressure fluctuation at a distance r away from the 
monopole can be expressed by 

Q(t - — ) 

P - P - Pq - ^ (2.1) 

Unr 

where p is the pressure perturbation (* p - P 0 ) • Notice that p is 
proportional to the value of Q at the retarded time (r = t - r/c Q ) , 
i.e. at the instant when a wave travelling at speed of sound c Q had to 
be launched in order to reach the observer at the current time t. 

The next simplest elementary solution is a dipole, where sound is 
generated by the injection of momentum. In other word, the acoustic 
dipole is equivalent to a concentrated point force, varying in 
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magnitude or direction or both. The strength of a dipole is the force 

itself. The pressure fluctuation seen by an observer 

at x - (x^,X2,X2) relative to the dipole can be expressed by 


P 



( 


F l< c 



4?rr 


) 


( 2 . 2 ) 


where r — |x| . It should be noted that the dipole strength, as a 
vector, has direction as well as magnitude. 

Gutin (1936) was the first to consider flow- generated acoustics 
in his pioneering study of propeller noise. He applied the results of 
Lamb (1932) on the acoustic field of a concentrated force (i.e. 
dipole) to the problem of sound radiation from a stationary propeller. 
Yet, it was not until Lighthill (1952) introduced the acoustic 
analogy, then a general theory of modern aeroacous tics began to 
emerge. Lighthill 's acoustic analogy will be derived and discussed in 
the following section. 


2.2 GOVERNING EQUATIONS FOR MODERN AEROACOUSTICS 

In this section, we will derive, review and discuss the 
generalized Lighthill 's equation and an extension to it by Ffowcs 
Williams and Hawkings (1969a). 


2.2.1 Generalized Lighthill f s Equation 

Modern aeroacoustics owes its major development to Lighthill 
(1952), who first considered turbulent flow as a source of sound. 
Originally Lighthill 's theory was introduced to calculate the sound 
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radiation from a relatively small turbulent flow region embedded in an 
unbounded homogeneous fluid at rest, where the speed of sound c Q and 
the density of the undisturbed fluid pQ are constants. Away from the 
small region of turbulence, the density fluctuation p (= p - p Q ) can 
be considered small and varying isentropically , hence it satisfies the 
homogeneous wave equation 

A- -4 A--* 

at ax t 

The pressure fluctuation p (- p - P Q ) can be expressed in terms of the 
density fluctuation p according to the isentropic relationship 

P - P 0 ~ °0 G ' V (24) 


Ligh thill ' s analogy can be derived from the equations of mass and 
momentum conservation, if they are arranged in a form similar to 
equation ( 2 . 3 ). We start the derivation of the Lighthill's equation 
by writing down the continuity equation (in Cartesian tensor notation) 


dp 8 

— - — + 

at dx L 




Q 


( 2 . 5 ) 


where Q is the rate of introduction of mass per unit volume. The 
momentum equation can also be written in the same fashion as 


A <?u i > * 4 


d °ii „ 

+ F. 

ax, 1 


( 2 . 6 ) 
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w here F ] _ is the external force per unit volume acting on the fluid, 
anc * a ii stress tensor given by 


a ij = - w u + 


( du i 3u j 2 5u k I 

^ i a Xj + dx. ' 3 5 ij d^ } 


where <? is the hydrostatic pressure (assuming Stokes' hypothesis 
holds). Eliminating pu^ between equations (2.5) and (2.6) gives 


dp 8 Q dF i d l 

8t 2 dt dx . dx dx 

1 i j 


(?u i u j - a ij) 


(2.7) 


Substracting d (Cq p 6^j)/dx^dXj from both sides of equation (2.7), 
Lighthill obtained an inhomogeneous wave equation which is later 


refered to as the Lighthill 's equation (or the Lighthill acoustic 
analogy) : 


2 2 2 
8 P 2 3 P <9Q 3F i 8 T. . 

~2 ~ C 0 2 + ~ 

dt dx. dt dx. dx dx 

1 1 i j 

where Th. is the Lighthill 's stress tensor defined by 


T ij - ?“i u j 


”ij 


2 

c p6 ■ ■ 
(T ij 


- Pu.u. 


2 , du. 

hj <p ’V ) • “ Inj 


du . 


3x i 


2 du k 

3 ( d^^ij ) 


( 2 . 8 ) 


(2.9) 
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The three terms which appear on the right hand side of equation 
(2.8) are generally recognized as monopole, dipole, and quadrupole, 
respectively . 

From equation (2.9), one can easily see that the quadrupole 
sources represented by Lighthill's stress tensor , actually consist 
of three different mechanisms: the first term pu-Uj represents the 

direct convection of the momentum component pu^ by the velocity 
component Uj . the second term fi^p - p) represents the effects of 

any change in entropy, and the last term simply represents the viscods 
stress effects (for an ideal fluid, this represents the convection of 
momentum by the motions of molecules relative to the flow velocity 

up. 

In Lighthill's original paper, was approximated by Pg u i u j ^ or 

turbulent flow at small Mach number assuming the density fluctuation p 
is much smaller than p Q ; under this condition, the flow can be 
considered to be isentropic (i.e. p - Cq p - 0) , thus the second term 
in equation (2.9) can be dropped. Also the viscous effects will be 
small compared to the inertial effects (assuming Re » 1), hence the 
third term in equation (2.9) can also be dropped. 

In addition, Lighthill assumed that neither mass sources nor 
external forces are found in the acoustic medium (thus Q - - 0) . 

Thus the original Lighthill equation was given by 

2 2 2 
dp 2 dp d T i j 

o " c 0 9 ” I 

dt l 3x^ 3xpx 


-24- 


Clearly, the Lighthill equation (2.8) is an exact equation. No 
assumption was made during its derivation from the continuity and 
momentum equations. Also it should be noted that the Lighthill 
equation applies only to fluid without internal solid boundary. The 
theory for sound generated aerodynamically with the presence of 
internal solid boundaries was later developed by Curie (1955) and 
Ffowcs Williams and Hawkings (1969a), this is the subject we will 
discuss in the following section. 

2-2.2 Theory for Flows with Internal Boundaries 

From the previous section, we had shown that the problem of 
aerodynamic sound can be posed as an acoustic analogy in which the 
turbulence provides a quadrupole distribution in an ideal acoustic 
medium at rest. However, Lighthill 's general theory did not consider 
the general effects of any internal solid surface that might be 
present. The first attempt to extend Lighthill 's general theory was 
done by Curie (1955). Curie considered two effects introduced by the 
solid boundaries: (1) the reflection and diffraction of sound wave at 
the solid boundaries, and (2) the pressure dipole distribution on 
solid boundaries (it is the limiting case of Lighthill 's quadrupole 
distributions). Curie showed that the effect of the solid boundaries 
upon the sound field is equivalent to a surface dipole distribution, 
each representing the force that solid surfaces applied to the fluid. 

In the hydrodynamic sense, a full representation of any solid 
surfaces lying within the fluid should consist of distributions of 
sources and surface forces. Thus, later Lighthill (1962) further 
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suggested that the sound field should be described by surface 
distributions of monopoles, dipoles, and distributions of quadrupoles 
in the space exterior to the solid boundaries. In this section, our 
review will closely follow the study of Ffowcs Williams and Hawkings 

(1969a) . 

It is quite clear that the theory of aerodynamic sound is built 
upon the equations of mass and momentum conservation of a compressible 
fluid, and results in an inhomogeneous wave equation governing the 
production and propagation of sound waves. These equations are 
generally valid only in the region exterior to any closed internal 
surfaces that may be presented. Such a situation is essentially 
inhomogeneous in space, in that these equations are valid in the 
volume outside the surfaces, but are meanless elsewhere. However, the 
spatial homogeneity can be restored by considering the following 
problem. First the closed solid surfaces are replaced with their 
corresponding mathematical surfaces, thus instead of having solid 
surfaces embedded in the fluid, the fluid is petitioned into several 
regions. The motion of the fluid on and outside the mathematical 
surfaces is defined to be completely identical with the real motion, 
whereas the interior flow can be specified arbitrarily. The interior 
motion is usually assumed to be very simple, and consequently does not 
match the exterior flow at the boundaries. Mass and momentum sources 
then have to be introduced to maintain these discontinuities, and 
these ultimately act as sound generators. 

Consider an unbounded fluid with a moving solid surface embedded. 
Replacing the moving solid surface by an equivalent mathematical 


surface 
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f(x,t) - 0 


( 2 . 10 ) 


f > 0 represents the volume exterior to the surface, and f < 0 denotes 
the region enclosed by the surface f - 0. Let v be the velocity of 
the solid surface f, and n be the outward normal on surface f — 0. 
Assume that in the region inside f - 0 the fluid is at rest (i.e. ^ - 
0) , with density Pq and pressure Pq . These values of density and 
pressure are those which would be found in the real fluid were it at 
rest. Also assumed is that the surface is impermeable, i.e. Ui - v. 
on the surface f - 0 . The continuity and momentum equations can then 
be written as (Ffowcs Williams and Hawkings, 1969a) 


dp 

dt 



P o v i 5(f) 



rt (?u i > + 


577 au i u j > 


axj 
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5(f) 


df 

axj 


(2.11) 


Similar to what we have done in the previous section, pu. is 
eliminated from the equations. Thus we obtain the Ffowcs Williams and 
Hawkings ' equation 


a 2 . 


at 2 


.2 

2 dp 


'0 2 

dx L 

i 


dx^dXj 


dx 




£f 

dx 


) + 


( 2 . 12 ) 


It should be noted that the Lighthill's stress tensor is zero for 

fluid inside surface f. 

Equation (2.12) shows that, in general, sound can be regarded as 
generated by three source distributions. The first is a distribution 
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of quadrupoles of strength , distributed throughout the region 
exterior to the solid surfaces (Lighthill, 1952). This is 
supplemented by surface distributions of dipoles of strength "ij n j 
(Curie, 1955). If the solid surfaces are moving, the sound field will 
be further enhanced by a distribution of monopoles of strength p Q v i n i 
to represent the volume displacement effects of moving solid surfaces 
(Ff owes Williams and Hawkings, 1969a). 

It should be noted that equation (2.12) remains valid even when 
there are shock discontinuities embedded in the flow field. The shock 
surfaces can be treated just like solid surfaces by replacing them 
with the corresponding mathematical surfaces; since the mass and 
momentum fluxes are continuous across shock, no distribution of 
monopole or dipole is needed on the shock surface. But the 
distribution of ■ will now contains discontinuities other than at 
physical boundaries. 

2.3 GENERAL SOLUTION TO THE GOVERNING EQUATIONS 

No matter what form of governing equation is chosen, we generally 
have a nonhomogeneous wave equation to solve. Consider a general 
expression for such type of partial differential equation 



2 

2 d p 



g(x, t) 


1 

where g(x,t) represents 
combination of monopoles 


the forcing terms which are 
dipoles, and quadrupoles. 


(2.13) 


the linear 
Refer to the 


right hand side of equation 


(2.8) for details . 
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The most straightforward way to solve a nonhomogeneous wave 
equation is by using the Green's function method. The Green function 
of wave equation is essentially the solution for a nonhomogeneous wave 
equation with an unit impulse as its forcing term. Let y = (y^.y^’Yg) 
be the source coordinates where an impulse of unit strength is applied 
at the instance when t “ r. The Green’s function of wave equation 
G(x,t;y,r) satisfies 



9 

2 d^G 
: 0 ? 

l 


S (x-y) 


S(t-r) 


(2.14) 


G must also satisfy the causality conditions for hyperbolic equations , 
i . e . 


G 


dG 

Jt 


0 ; 


for t < r 


(2.15) 


G can be solved by first finding the Laplace transformation of 
(2.14) with respect to time, solving the resultant equation by use of 
the properties of spherical symmetry, and then inverting the 
transformation to find the solution (Carrier and Pearson, 1976) 

&(t - t + r/c Q ) 

G 2 (2.16) 

47rc r 
0 

where r - |x-y| • The formal solution of (2.13) can then be expressed 
in the following integral form 



-29- 


+co 


5(r - t + r/c n ) 
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V -00 


g * 


dr 


(2.17) 


Equation (2.17) can be integrated by using the property of the 
Dirac delta function 


4-oo 

| f(x,t) S(t-r) dt - f(x,r) ( 2 ■ 

- 00 

By first carrying out the integration in (2.17) over time (r) using 
(2.18), we can have 



where the square bracket indicates that g and r are evaluated at the 
retarded time r (■ t - r/c Q ) for each point in the volumetric 
integration . 

Replacing g in equation (2.19) with the right-hand side of the 
Ffowcs Williams and Hawkings’ equation (2.12), its solution can be 
found as 
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( 2 . 20 ) 
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v/here v r is the normal component of the velocity of the solid surfaces 
S, V represents the volume exterior to S, and M r is the component of M 
in the radlatlon direction. Note that the terms in the right-hand 
side of (2.20) are integrated with respect to the Lagrangian 
coordinates rj which move with the acoustic sources. The second and 
the third terms on the right-hand side of (2.20) are linear, they are 
generally well defined and comparably easy to evaluate (surface 
integrals). However, the proper treatment of the quadrupole term 
requires integration over a volume that is large enough to enclose all 
nonlinear effects. The computation will clearly be very lengthy. 
Fortunately in most low speed cases, quadrupoles are not as efficient 
sources as monopoles and dipoles, i.e. the quadrupole terms are often 
negligible . 


2 SOUND FIELD OF MOVING SINGULARITIES 

According to Lighthill (1962), the sound field for a point force 
in uniform rectilinear motion can be given by 
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z i "Yi 


3F 
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3 2 2 

L 4jrc Q r (1-M r ) 3t 


( 2 . 21 ) 


where and y^ are the Cartesian coordinates of observer and source 
respectively. This equation implies that if the force is constant no 
sound will be radiated. However, it is generally known that sound can 
also be generated by the convection of constant forces. The uniform 
rectilinear motion assumption made in equation (2.21) is sufficiently 
restrictive to remove some of the terms giving rise to the radiated 
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sound. This effect has been studied by Lowson (1965); in this section 
our discussion will closely follow his derivation. 


2.4.1 Sound Field for a Point Force in Motion 

Consider a point force in arbitrary motion. It can be expressed 
in terms of the product of F and a Dirac delta function 5(y). F 
gives magnitude as well as direction of the point force, and it can be 
regarded as a function of time only. 8 - 6(y) gives the position of 
the moving point force. Since the point force is moving, y is also a 
function of time. 

By putting g - -3F i /dx i into equation (2.19), the sound field of 
a moving point force (dipole) can be written as 


P 



. S(F . 8 ) 
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r 3x^ 
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( 2 . 22 ) 


The proper evaluation of (2.22) depends critically on the 
the square brackets, which requires the evaluation of its 
the retarded time r (= t - r/Cg). 

Let us consider any function f(x,t). The chain Rule 


treatment of 
contents at 

shows that 
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For a moving singularity, we have 
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(2.24) 
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where y i is the component of y in the i direction. From (2.23) and 
(2.24), with some arrangements, we have 


r ds “ 

i r c o M i ' 



L at _ 

[ J 

2 L o j 

ax i 

(2.25) 


where M r is the component of the instantaneous convection Mach number 
in the direction r toward the observer. M r car. be written as 

( 2 i -y i ) Mi 

M r - — r 


is the component of the instantaneous convection Mach number 
in the i direction. 

Using (2.23) and (2.25), after some mathematical manipulations , 
equation (2.22) can be rewritten as 



(2.26) 


where [5] - 5 (z-y(t-r/c Q ) ) . To evaluate the integral of [5], the 
variable in the integrand is changed from y to r) 


n - z - y(t - r/c Q ) 


Thus the integration can be evaluated directly, yielding 
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Thus, the radiative terms of (2.28) can be given as (Lowson, 1965): 


(2.28) 
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where aM r /at is the component of acceleration in the r direction. 

3M r (z i -y i ) 

It r at 

The second term in (2.29) will vanish if the convection velocity 
is constant, and the result reduces to the same as given by Lighthill, 
i.e. equation (2.21). It should be noted that if the force is 
accelerating, the second term will give rise to the sound radiation 
even when the force is constant. We can also write down the near 


field term, which can be given as (Lowson, 1965). 
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Note that the near field solution shows no dependence on the 
acceleration. 


2-4.2 Sound Fields for Moving Monopoles 

The classical solution for a monopole with constant convection 
velocity can be written as (Lighthill, 1962) 


1 3Q ' 
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Just as we have done for the moving dipole, the far field 
solution for a moving monopole can be written as (Lowson, 1965) 
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The near field solution can also be found as (Lowson, 1965) 


Q(M r -M 2 ) 

. 4*c 2 (l-M r ) 3 r 2 . 


(2.33) 


For a monopole with constant convection velocity, the second term 
in (2.32) vanishes. However, this is different from the classical 
result (2.31) by a factor of (1-M r ). To clarify this, we rewrite 


equation (2.32) as 
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We can now see that the sound field of an arbitrarily moving 
monopole is actually caused by three different mechanisms: (1) pure 

monopole effect of the rate of introduction of mass, this yields the 
same result as the classical solution, (2) dipole effect of the 
convected momentum, and (3) quadrupole effect of acceleration. 

2.5 SOUND FIETJt OF SINGULARITIES IN CIRCULAR MOTION 

The theory of the far field solution for rotating singularities 
was studied by Ffowcs Williams and Hawkings (1969b). In this section, 
their results will be reviewed. Unlike the method we discussed in 
section 2.4, which deals with the solution in time domain, Ffowcs 
Williams and Hawkings solve the problem in frequency domain. 

From the discussions in previous section, the solution of (2.13) 
can be written in real time form as 

p (x , t) f g(y.t-r/c 0 ) dy (2.35) 

4 * c 0 r V 

alternatively it can be written in spectral form as 


- 2irif r/c 
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g(- - r.f) 
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4ttc r 
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p(x.f) “ 


(2.36) 
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where r is the unit vector in the radiation direction. The solution 
in spectral form p(x,f) is the generalized Fourier transform of p(x,t) 
and can be written as 


+00 


p(x, f) 


I 


p(x,t) e ' 27rlft dt 


g(k»f) is the wave number, frequency spectral component of the 
source field, it is the four-dimensional Fourier transform of g(x,t) 


gO, f) 


-H=o 


If 


g(x,t) e 


-2xik*x -2xift 


dxdt 


(2.37) 


In terms of the spectral components of the multipole strength of 
equation (2.8), this can be written as 


g(k,f) - S (k , f ) + 2xik i D i (k,f) 


(2x) 2 k i k j T ij (k,f) 


The particular component of g(k,f) that generates acoustic waves 
at frequency f, propagating in the direction r, is 


g(- — r,f) 


S(- — r,f) - 2xi — r, D.(- — r,f) 
c 0 c 0 1 1 c 0 


- (2*) 2 ^ T l (- r , f ) 

c^ J c 0 


(2.38) 


where ^ is the component of the unit vector r in the direction i. 



-37- 



Figure 2.1: Sketch of the Rotating Acoustic Singularities. 

Consider concentrated acoustic sources rotating in the *3 " 0 
plane at radius R with rotational frequency fl. The acoustic sources 
can be written symbolically as 

g(x,t) - , t) •& (x-R) 

where q(<£,t) is the time dependent strength of source which occupied 
angular position <f> from the x axis at time t ~ 0, and R is the 
position vector of that source at time t. Refer to Figure 2.1 for the 
geometry. Putting this into equation (2.37), the four-dimensional 
Fourier transform of the source density becomes 
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+<o 

g(k.f) - J q(*,t) e - 2,rik,R e ' 2nift dt 

- qo 

Again q(0,t) may be expressed in terms of its generalized Fourier 
transform as 


+oo 

q(^.t) - J q(^,a) e 2wiat ^ 

- 00 

and the radiating component of the source spectrum becomes 

+co 

g(- ~ r,f) - JJ q(<A ( a)e- 27 rit < f - a )e 2,rifRsintfcos ( 2,Tnt+ ^)dtda ( 2 . 39 ) 

^ - oo 

In (2.39), r.R is equal to Rsin0cos(2*nt+^) , and 8 is the angle at 
which the wave is propagating relative to the rotor axis x 3 (which 
is 8 - 0°) . 

To evaluate equation (2.39), the integration with respect to time 
is performed first, yielding 


+«o n-+<x> 


g(- — r , f ) - f Y q(*, Q ) e -W2 e i^(f-Q)/0 J , 2xfRsintf. 
c 0 J i—. n v c 


-) 


- oo n- - 
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Y q(rf,f-nO) 

n--® 


e - in(n/2-<f>) j ^ 2wfRsin 8^ 


(2.40) 



-39- 


where J is the Bessel function of the first kind (subscript n denotes 
n 

the order of the Bessel function). The infinite summation over the 
Bessel function is due to the Doppler effect of rotating 
s ingularities . 

Consider the three different source components: monopoles, 
dipoles, and quadrupoles . Their spectra can be written as s 
d i (^,a), and tjj (<p ,a) respectively. By equations (2.40) and (2.36), 
the far field general solution for p in spectral form can be expressed 


as 


- 27 rifr/c 


p(x, f) 




n~ r 

- !{(•«•*• 


nil) - r . d. W , f -nfl) 


o n— « 


- (2*) 2 rj.rjtijW.f-nn) ) 

°0 

-in(w/2-*) j (. 2fffRsin ^- ) 
n c 0 ^ 


(2.41) 


To make equation (2.41) practically useful, one must find the 
power spectral density relations to predict statistics of the 
radiation field in terms of statistics of the sources' strength. This 
is done by multiplying both sides of equation (2.41) by their 
conjugates, averaging, and normalizing by a factor common to both 
sides whose .magnitude need not be considered and can be taken as 
unity. Since our present study is dealing primarily with dipole 
radiation, we limit the discussion to this particular case. Equation 


(2.41) reduces to 
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(2.42) 


This equation will be used to develop a formula relating the power 
spectral density of a time -dependent dipole strength to the power 

spectral density of the density fluctuation in the far field radiation 
field, 9(x,f) 
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(2.43) 


Here d^, has been written for the dipole component in the 
particular radiation direction r^d^. If we assume that all Fourier 
elements are uncorrelated unless they are conjugates, and this is so 
in a stationary field, we can have 


d r (^,a)d r (<f> 9 0) - D r (^,a); if a - 


- 0; 


if a / -/? 


(2.44) 


Furthermore, we can neglect the dependence of D r on a without 
losing any generality, since 0(x,f) is the averaged quantity over one 
revolution period. Thus the double sum in equation (2.43) can be 
reduced to a single summation over n: 
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where is the power spectral density of the dipole strength in the 
radiation direction. 

The spectra for pressure fluctuation p(x,f) can be obtained from 
the spectral for density fluctuation p(x, f) in conjunction with the 
isentropic relation given in equation (2.4): 

p(x,f) - c^ p(x,f) 


Thus, the power spectral density of the pressure fluctuation 
<S(x,f)> can also be found as 

<S(x,f)> - 0(*.f) / p Q c Q 
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(2.46) 


This is the primary result of this section. It is extremely 
useful in the prediction of high frequency broadband rotor noise. The 
only term that remains unknown is D r> the power spectral density of 
the dipole strength in the radiation direction. It can be found from 
the unsteady airfoil response or surface pressure fluctuations due to 
various blade - turbulence interactions and will be discussed in detail 


in Chapter III . 
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2 • 6 APPROACHES TO THE PREDICTION OF BROADBAND ROTOR NOISE 

In analyzing the broadband rotor noise problems, several 

approaches have been applied so far by previous researchers. They 
are : 

(1) Treat the general case of unsteady forces distributed in 
space following equation (2.8) with specialization to 
rotating blades . 

(2) Approximate the distributed blade unsteady forces as rotating 
point forces (i.e. dipoles) using the method we just 
discussed in section 2.5 (Ffowcs Williams and Hawkings, 

1969b) . 

(3) Find the acoustic radiation from an airfoil segment in linear 
motion, and then approximate a rotating blade section by a 
sequence of airfoil straight line motions. 

An example of the first type of approach is the analysis of 
Homicz and George (1974). The second type of approach has been used 
by George and co-workers on the study of a number of broadband rotor 
noise mechanisms (George and Kim, 1976; Kim and George, 1980; George, 
Najjar, and Kim, 1980; Chou and George, 1985). Examples on the third 
approach are the analyses of Amiet (Amiet, 1976; Schlinker and Amiet, 
1981). A brief review of some of the representative analyses is given 


below . 
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2.6.1 The General Approach 

The analysis of the sound generated by the unsteady forces due to 
^urbulewce injection into subsonic unducted rotors began with the 
quite general analysis of Homicz and George (1974). In their 
analysis, they found that the far field sound from arbitrarily- varying 
forces on a plane S (e.g. a rotor disk) can be expressed as 


<S(x,f)> 


x.x. 


4 pc* 

0 0 


91 


<P L L ^ ’ 

i j 


f "o^ dS 
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(2.47) 


where <P T (x,k,f)> is the generalized averaged three-dimensional 

l j 

cross-spectral density of the forces (e.g. rotor loads) per unit 
area on the plane S (assume to be x^ - 0) , and tiq is the unit vector 
along the observer's direction. Refer to Figure 2.2 for the geometry 
of the analysis. 

It should be emphasized that P L>L is the load spectrum as 
measured in coordinates fixed with respect to the observer, not moving 
with the blades. Thus blade-to-blade correlation can be treated 
rigorously . 


2.6.2 The Rotating Dipole Approach 

The analysis of George and Kim (1977) was developed for the high 
frequency inflow turbulence noise prediction. Their analysis is based 
on the second approach mentioned previously. 

By using the rotating dipole theory of Ffowcs Williams and 
Hawkings (1969b), they assumed the fluctuating forces rotate in a 
circle and their components in the observer's direction are 
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v. Observer 



Figure 2.2: Plane Rotor Geometry and Coordinates System. 

statistically stationary. This effectively restricts the analysis to 
forces that are normal to the rotor plane (fluctuating lifts) and thus 
does not account for the smaller fluctuating drag forces which can be 
significant for observers in the rotor plane. Also, since only dipole 
forces are assumed, the analysis does not account for the detailed 
radiation directivity of the blade elements. The assumptions of 
circular motion also does not allow accurate treatment of the 
helicopter forward- flight cases. Fortunately, these restrictions are 
not very important except for observers within about 15 degrees to the 
rotor plane or for rotor advance ratios greater than about 0.4. 

This approach was later used by George and co-workers on other 
rotor broadband mechanisms including boundary layer/trailing edge 



-45- 


noise (Kira and George, 1980), tip vortex formation noise (George, 

Najjar, and Kim, 1980), and turbulent vortex shedding noise from blunt 
trailing edges (Chou and George, 1985). This approach generally gives 
a reasonably simple equation and can be evaluated quite easily, even 
at high frequencies* 

2.6.3 The Computational Approach 

The rotor noise prediction methods of Amiet (Amiet, 1976, 

Schl inker and Amiet, 1981) can be classified as the third approach we 
mentioned in the beginning of this section. His analysis is a 
semi-analytical, semi-computational approach. The circular motion of 
a rotor section is approximated by a sequence of an airfoil segment s 
linear motions. The noise spectrum of an operating rotor can thus be 
obtained by summing and averaging the noise spectrum of a number of 
individual blade linear motions over one revolution. A schematic 
sketch of this approach is given in Figure 2.3. 

Amiet' s method has the advantage of being able to account for the 
full range of wavelength-to-chord ratios and thus accurately predict 
the directivity of the sound radiation. In addition, it is very easy 
to account for the effect of forward flight. The equation is simple 

and is very efficient to compute. 

In this chapter, we have reviewed the analytical tool for 
calculating the far field noise from given sources. The rest of this 
thesis will deal with the physical modelling of various rotor 


aerodynamic noise sources. 



Chapter III 

ROTOR BROADBAND NOISE MECHANISMS 

In the present chapter, various "first-principles" analyses of 
rotor broadband noise are reviewed and extended; new analyses are 
developed. Emphasis is placed on the physical modelling of various 
broadband noise mechanisms. The results are used in association with 
the basic acoustic analysis given in chapter II to predict the 
absolute noise spectra. These analyses are not based on empirical 
correlation equations and do not require the determination of 
empirical constants for different families of rotors. To determine 
the relative importance of various mechanisms for different rotor 
parameters and for different operating conditions, computations are 
made for many rotors. The results are then compared to each other and 
to selected experimental data. 

3.1 REVIEW OF PREVIOUS RESEARCH 

Broadband rotor noise had been studied extensively for almost two 
decades. Historically, the first broadband noise prediction methods 
were based on empirical correlation of overall sound pressure level 
(OASPL) , e.g. Widnall (1969). Very early investigators identified 
broadband noise with some sort of vortex shedding from rotor blades 
(similar to that from circular cylinders), hence the early name 
"vortex noise" was used. Actually most broadband noise is due to 
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unsteady loading and surface pressure fluctuations induced by 
turbulent flows . 

Vojetex shedding remains as one noise source; it was found that in 
the atypical case of laminar flow, the laminar boundary layers on 
blades can indeed shed nearly regular vortices from trailing edges and 
thus radiate a narrow peaked broadband sound, this was refered to as 
"high frequency broadband noise" by several investigators (Paterson et 
al., 1973, Aravamudan et al . , 1979a). This mechanism is not important 
for most full scale rotors since the flows around them are generally 
turbulent. Possible exceptions are helicopter tail rotors and small 
axial flow fans. Even in these cases it can be eliminated easily by 
tripping the boundary layers (Aravamudan et al., 1979a, 1979b), 

Reviews of recent broadband rotor noise research can be found in 
the work by Schlinker and Brooks (1982), or by George and Chou (1983a, 
1984a) . 

3-1.1 Inflow Turbulence Noise 

Inflow turbulence noise originates from unsteady loading 
(particularly lift) fluctuations when rotor blades encounter 
turbulence. The noise produced can be broadband, or, if turbulent 
eddies are chopped by more than one blade, the noise can be 
narrow-bank random, i.e. spectrum humped or peaked around harmonics of 
the blade passage frequency. 
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3 . 1 . 1.1 Analysis of Honicz and George 

The analysis of the sound generated by unsteady loading 
fluctuations due to turbulence ingestion began with the quite general 
analysis of Homicz and George (1974). They treated the general case 
of unsteady forces distributed in space, following the Lighthill 
equation of aeroacoustics with specialization to rotating blades. 
Refer to Section 2.6 for the acoustic result of their analysis. 

Homicz and George then applied this acoustic result, i.e. 
equation (2.47), to the unsteady loading (lift) fluctuations on 
rotating blade. The unsteady loadings were obtained from an 
approximate compressible aerodynamic analysis for an inflow of 
homogeneous isotropic turbulence. The noise spectrum they obtained 
can be written as 
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E, 


„ 2 , t 2 


int “ sin * J nB-l (M 0 a cos ^ ) 


(^)V-(l -) 2 


"turb 


(1 + 4 * 2 £ 2 ) 3 


where Mq, M c , and M fc are the rotational, convective, and turbulence 
Mach numbers, and £ - kA is the dimensionless turbulence wave number. 
R 0 represents the effective radius of the equivalent rotating point 

dipole. The notation nB denotes the nearest blade passing harmonics 
to f/Q. 

It should be noted that while equation (3.1) was derived for 
inflow turbulence defined by Dryden spectrum, other homogeneous and 
isotropic spectra can also be handled, since the integration over £ 
must be done numerically in any event. (In such cases, E Curb will have 
to be modified . ) 

To account for the effect of compressibility in the unsteady 
aerodynamic analysis, Homicz and George found that an additional 
multiplicative factor should be added to the integrand in equation 
(3.1), which is 


E 


aero 



(3.2) 


where 
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E t . (Rj ♦ E|)V2 


and Ic and tc are the averaged reduced aerodynamic frequencies along 
c s 

the chord and the span directions over one blade revolution. 

Equation (3.1) was derived assuming uniform spanwise source 
strength. This is obviously not the case if the turbulent wavelength 
is smaller than the blade span. In order to account for this effect, 
Homicz and George (1974) suggested that another multiplicative factor 
should be added to the integrand in equation (3.1): 


E 


span 


(1 + E c b/frc) 


(3.3) 


It should be pointed out that the analysis of Homicz and George 
uses coordinates fixed with respect to space rather than to the rotor 
blade. This allows rigorous treatments of the blade-to-blade 
correlation. Such correlation explains the humped or peaked nature of 
the low frequency part of the spectrum as being due to the large scale 
components of the turbulence inflow. These large scale components 
give nearly periodic (i.e. modulated) disturbances as they are swept 
through the rotor plane; this leads to a nearly periodic but finite 
bandwidth sound. The most important parameter affecting the shape and 
magnitude of the acoustic spectrum is shown to be the ratio of the 
time taken by the rotor to complete one revolution (1/fl) to the time 
needed to convect one integral length scale of turbulence through the 

rotor disk (A/V c ) . 
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Thls analysis is not well suited for high frequencies, where the 
complexity of combining equations (3.1), (3.2), and (3.3) requires 
prohibi-tly large CPU times for the calculations. Thus, high frequency 
analyses were developed later by George and Kim (1976) and by Amiet 
(1976). 


3. 1.1. 2 High Frequency Analysis of George and KLm 

In order to predict the high frequency noise efficiently, George 
and Kim (1976) treated the same problem in a different fashion. They 
approximated the distributed blade forces as rotating concentrated 
forces (dipoles) , and used the result of Ffowcs Williams and Hawkings 
(1969a) to get the noise radiation. The result of their analysis can 
be written as 
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It should be noted that <j> is the observer elevation angle 
measured from the rotor plane, as shown in Figure 2.2. Equation (3.4) 
completely neglects any effect on blade-to-blade correlation, and is 
only valid at high frequencies. This analysis also assumes that the 
force components in the observer's direction are statistically 
stationary. This effectively restricts the analysis to the forces 
normal to the rotor plane (i.e. lift) and thus does not allow 
accounting for detailed radiation directivity nor for the much smaller 
torque forces which can be significant for observers near the plane of 
rotation. The assumption of circular motion also does not allow 
accurate treatment of forward- flight helicopter cases. However, these 
are not important restrictions except for within about 15 of the 
rotor plane, or for advance ratios greater than about 0.4 (George and 
Chou, 1983a, 1984a). 

This formulation can be evaluated quite efficiently since the 
integration inside the summation is evaluated semi-analytically . 
Details on the numerical approach are given by George and Kim (1976) . 

The original analysis was specially derived for turbulence 
defined by a Dryden spectrum. Other homogeneous, isotropic turbulence 
models can also be handled by modifying , the power spectral 
density of the lift on rotor blade. But, the equation is then more 
difficult to evaluate compared to the original equation, which can be 
evaluated quite efficiently. 
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3. 1.1. 3 Analytical-Computational Analysis of Amiet 

The analysis of Amiet (1976) is based, on a different concept. 
Initially, Amiet analyzed the radiation of sound from a stationary, 
non- rotating airfoil in a uniform mean flow containing turbulence 
(Amiet, 1975; Paterson and Amiet, 1977). Amiet’s analysis can handle 
a full range of wavelength-to-chord ratios, and can accurately predict 
the directivity of the sound radiation (which becomes multi-lobed at 
intermediate wavelength-to-chord ratios). Figure 3.1 shows several 
representative directivity patterns (Kaji, 1975). The result is given 
by Amiet as 
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Refer to Amiet' s original paper (1975) for details of this analysis. 

Later, Amiet used this result to synthesize the average noise 
radiation from the rotating blade by numerically summing and averaging 
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Figure 3.1: Comparison of Directivity Patterns of Radiated Sound 

Predicted by Compact and Non-Compact Sources, M - 0.5 
from Kaj i (1975) . 


the radiation from a sequence of blade linear motions which 
approximate the circular (or epicycloidal) motion of a rotating blade. 
Refer to Figure 2.3 for sketch of Amiet's approach. 

The analysis of Amiet includes both low- and high-frequency 
response functions in addition to source non -compactness and 
compressibility effects. Amiet’s method has the advantage of being 
able to treat forward- flight easily and of being based on a more exact 
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model of sound radiation directivity. However, when one sums and 
averages the multi-lobed radiation pattern over one revolution of the 
rotor blades, the pattern is smoothed out to a pattern which, except 
when the observer is very close to the rotor plane, is quite close to 
the results of the approximation of George and Kim (1976). 

3.1.2 In- Plane Broadband Noise Due to Turbulence 

The analyses described in the previous section considered lift 
fluctuations due to turbulence ingestion as the primary source. They 
tend to under-predict the noise perceived by an observer in the plane 
of the rotor, where the smaller drag dipoles and other in-plane 
sources become dominant. 

Consider an in-plane fluctuating force D (e.g. drag). According 
to elementary acoustics, the sound field can be written as 

. 1 d f f °/9 

4* dx 0 J [ r(l-M r ) J (3.6) 

where 0 - 1,2 are the two in-plane directions and £ is the mean blade 
planform . Unfortunately, it is very difficult to obtain the correct 
distribution for the high frequency drag force component, even if the 
viscous effects are neglected. This distribution for the high 
frequency drag force component depends directly upon the full pressure 
distribution over the blade surface. The finite- thickness and the 
three-dimensional effects must be retained in calculating such 
detailed pressure distribution. In other word, simple aerodynamic 
response theory will not give the information we need for evaluating 
the in-plane broadband noise due to fluctuating drag dipoles. 
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Consider a thin airfoil embedded in turbulent flow. In order to 
satisfy the unsteady boundary condition (including velocity 
perturbations due to turbulence) at the airfoil surface, an additional 
monopole distribution should be added. Hawkings (1978) showed that 
for a thin symmetric airfoil in incompressible flow, the sound 
generated due to the in-plane components of the velocity perturbation 
Wp can be given as 


p(x, t) 


1 d 2 
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Un dtdx^ 
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where /9 is again confined to the two in-plane directions . This 
mechanism is refered to as "unsteady thickness noise", and it has been 
shown to be a major in-plane rotor broadband noise source (Glegg, 

1986) . Note that this is the direct analog of thickness noise with 
the blade speed U replaced by the unsteady gust velocity w^. The main 
feature of this equation, however, is that it essentially represents 
an in-plane quadrupole field (because of the double derivative outside 
the integral) . 

Unsteady thickness noise is different from the noise generated by 
fluctuating drag dipoles. However, a simple dimensional analysis 
shows that unsteady thickness noise can be significantly higher than 
drag noise at high frequencies. A drag dipole typically has a local 
source strength of the order of O.lpUw (U is the typical blade speed), 
whereas the quadrupole source is of order of 0 . lpcw (c is the blade s 
chord). However, this has to be multiplied by the frequency f, 
because of the additional time derivative outside the integral. 
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Hence, the ratio of the unsteady thickness noise to the drag noise is 
of the order of fc/U . Clearly at high frequencies, the unsteady 
thickness noise will be much more important. 

Hawkings unsteady thickness noise analysis can be applied to the 
rotor broadband noise problem with an approach similar to Amiet's 
inflow turbulence noise analysis (1976). The sound field from an 
airfoil in linear motion must first be determined if this approach is 
to be used. 

Consider a symmetric airfoil whose surface is defined by 


y ~ T r 


{ »0<h 1/2 ♦ .pi) ♦ a 2 (i) 2 + a 3 (i) 3 ♦ . 4 (i) 4 ] 


Then the unsteady thickness noise can be given as (Glegg, 1986) 
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This analysis was derived with inflow turbulence defined by 
Karman spectrum. To find the noise due to this source, the sound 
spectrum is summed and averaged over one blade revolution. 

3.1.3 Boundary Layer/Trailing Edge Noise 

Noise can also be produced by the self -generated turbulence in a 
blade's boundary layer passing its trailing edge. This phenomenon was 
recognized as far back as 1959 (Powell, 1959). Various investigators 
looked at very simplified models for this noise, but these early 
models were not complete and were useful only as bases for empirical 
correlations. Fink (1978), for example, used such a correlation to 
predict the on-axis noise of a rotor due to boundary layer/trailing 
edge interaction. Complete "first-principles" analyses of rotor 
trailing edge noise were developed more recently by Kim and George 
(1980) and by Schlinker and Amiet (1981). 

The analytical problem of sound radiating from the effect of 
turbulence being convected past a non- rotating blade's trailing edge 
has been studied intensively since about 1970. A variety of model 
problems were studied (see the review of Howe, 1978), but these 
studies primarily resulted in scaling laws which require empirical 
constants be determined. There remained a number of questions 
regarding the details of the modelling and the effects of the Kutta 
condition. Alternatively, Amiet developed a method which is based on 
solving the problem of a statistically stationary pressure field being 
convected past a trailing edge (Amiet, 1976, 1978). Amiet's analysis 
requires that only the surface pressure spectrum in the boundary 
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layers be determined experimentally. Amiet's analysis has been 

compared to the experimental findings of Brooks (1981) and found to be 
consistent . 


3. 1.3.1 Analysis of Kin and George 

In 1980, Kim and George constructed an analysis of boundary layer 
noise from rotors by using the blade forces from Amiet's aerodynamic 
model in the same manner as they did earlier for the inflow turbulence 
noise. Thus, their analysis is again restricted to angles not too 
close to the rotor plane and to the low advance ratios which are found 
in helicopter forward flight. Again this is a high-frequency 
analysis, and it neglects blade-to-blade correlation. Their final 
result can be written as 


<S 1 (*,f)> 


Bf 2 bVsin 2 tf F g (|f-nO|)S pp (|f.nO|) 


o 3 2 
2np cr 

0 0 

j 2 (H 0 | cos*) 
n u ft 


I 


n ’-“ <1+ i 2 (|f-nO|) 


) (f-nfl)‘ 


(3.9) 


where 


F - F 2 + G 2 

O 


M+Mm+K^ 1/2 



1 

J [ (c 1 +s 1 )cos2K 1 +(c 1 -s 1 )sin2K 1 ) ]+l- (c 2 +s 2 ) 




A 

J [ (c 1 -s 1 )cos2K 1 -(c 1 +s 1 )sin2K 1 )]-(c 2 -s 2 ) 


c 1 - isj_ - E (2 /j( 1+M) 
c 2 - is 2 - E*(2(^+/iM+K 1 )) 
K l - wc/2U c 
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M - MK//9 2 

i 2 ( f) “ 2.1U c /2*f 

and the surface pressure spectrum S pp can be given as (assuming flat 
plate boundary layer) 


2 5* 2 x 10' 5 

pp u (1 + S + 0.2175 2 + 0.005625 4 ) 

u " 2*f 5*/U 
5 * _ 0.047 c Re~ 1/5 


, 0.1 < w < 10 


3. 1.3. 2 Analysis of Schlinker and Amiet 

Similar to his treatment to the inflow turbulence noise, Amiet 
first analyzed the sound field of turbulence past an airfoil's 
trailing edge (Amiet, 1976, 1978). He found that the sound field due 
to self- generated turbulence past the trailing edge of an airfoil in 
linear motion can be written as 

k Mcz 2 i 

<S 1 (x,f)> - [ T ) ^ |*(x,k x> A,k y ) | 2 i 2 (f) s pp ( f ) (3.10) 

v 47 rcr 


where S and are defined in the previous section. Refer to 

pp 2 

Amiet' s paper (1976) for details of this analysis. 

Schlinker and Amiet (1981) used the same numerical summing and 
averaging method that Amiet used for the inflow turbulence noise 
(1976) to treat the rotor trailing edge noise problem. Again, the 
simple rotating dipole method of Kim and George (1980) gives 
essentially the same results as Amiet 's method except within about 15 
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of the rotor plane, where additional source terms should be included 
in both methods. 

3-1.4 Tip-Vortex Formation Noise 

Another source of broadband noise on airfoils or rotors is that 
due to tip vortex formation. Similar to the boundary layer/trailing 
edge noise discussed in the previous section, it is generated due to 
surface pressure fluctuations beneath a highly turbulent tip vortex. 
Kendall (1978), Ahtye et al. (1979), and Fink and Bailey (1980) 
experimentally observed localized noise sources at wing and flap tips. 
Changes in noise from variations in rotor tip shape were 
experimentally observed some time ago by Lowson et al . (1972) 
although these effects may have been due to blade loading changes. 
George et al . (1980) identified this effect with the turbulence in the 
vortex formation and local separation region over the blade tip 
interacting with the trailing edge. 

The model starts with the experimental observations of separation 
on the suction side of blade or rotor tips due to the boundary layer 
being swept around the tip by the pressure gradient. This separated 
vortex flow is similar to the flow over the top surfaces of a 
sharp-edged delta wing in subsonic flow. 

Figures 3.2, 3.3, and 1.6 show the sketches of delta wing and 
wing/rotor tip vortex flows. It is generally known that these leading 
edge vortices are quite turbulent. Large fluctuating pressures have 
been measured on the surfaces of delta wings underneath these 
vortices. George et al. (1980, 1983b) used the pressure fluctuation 
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data from 2-dimensional flows, and data on the geometry and velocities 
associated with wing and rotor tip flows to estimate the separated 
turbulent pressure spectrum being convected past the trailing edge. 

This surface pressure information was used to predict the 
resulting radiated sound in a manner similar to George and Kim's 
treatment of inflow turbulence and boundary layer trailing edge noise. 
This tip vortex noise is shown to increase with blade loading, as had 
been experimentally observed in many cases. The original analysis 
gives the far field sound spectrum as 


<S, (*, f )> 


IlAVsin 2 * F g (|f-nO|)S pp (|f.nO|) 


« 3 2 

8irp c r 

0 0 


I 


n— co (1+ i 2 (|f. n n|) )(f ' nfl > 2 


J n (M 0 I cos *> 


(3.11) 


where 


s pp (f) - <0.5„ 0 v2) 2 L S l( 5> / V m 


u> - fL/V 
' m 


V m max iniura circumferential velocity in the vortex core, L is 

the size of the tip vortex separation region. S^S) is the normalized 
spectrum for the fluctuating surface pressure, which can be 
approximated by 
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c /->. n . for u> < 0.1375 

S^(cj) = 0, 

5 . 9703x10 “^c3- 3 . 5673x10 for 0.1375 < ^ < 0.387 

3 . 144 x 10 ' 3 sin(3.2388w-0.5506); for 0.3872 < 5 < 0.793 

1/(93. 035+557. 09w); for 0.7935 < w < 1.060 

l/( - 258 . 896+1964 . 19w-2416 . 78w 2 +1288 . 94w 3 -100 . 862w 4 ) ; 

for 1.0605 < cj 

Note that equation (3.11) is essentially the same formulation as 
the trailing edge noise analysis of Kim and George (1980) ; the only 
differences being the surface pressure spectrum used and a slightly 
different coefficient. (Also this mechanism occurs only at the tip 
region of a rotor blade) . 


3 2 INFLOW TURBULENCE NOISE 

It can be seen from Section 3.1 that the accuracies of the 
predictions of inflow turbulence noise depend upon how the inflow 
turbulence is described. To describe the inflow turbulence properly, 
one has to specify the turbulent intensity, the length scale, and the 

spectral shape. 

There are several equations available for describing turbulent 
spectra, e.g. Dryden spectrum, von Karman spectrum, and mild knee 
spectrum, etc. As pointed out by Houbolt (1973), nonhomogeneity can 
also change the spectral shape. Consider the time history shown in 
Figure 3.4. If the signals in two periods (a and b) were analyzed 
separately, spectral curves a and b would result, each of which would 
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be represented well by the von Karman spectrum. Analysis of the 
combined record, however, leads to the composite spectrum shown. This 
composure spectrum exhibits a "mild knee" variation, as contrasted to 
the sharper distinct knee of the Karman spectrum. 







Figure 3.4: Effect of Non-Homogeneity on Turbulent Spectral Shape, 

from Houbolt (1973). 

Also consider the data shown in Figure 3.5. To fit the same 
data, the von Karman spectrum yields a length scale of 700 ft, while 
the mild knee expression indicates a scale of 2100 ft. When one 
applies the two different spectra to the inflow turbulence noise 
calculation, the results can be quite different. 

In the present study, we will limit out discussion to the Dryden 
spectrum and the von Karman spectrum only. 
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Figure 3.5: Influence of Spectral Equations on Deduced A Values, from 

Houbolt (1973) . 


3 2.1 Dryden Spectrum and Karma n Spectrun 

The Dryden spectrum has been used by George and co-workers in 
inflow turbulence noise calculation (Homicz and George, 1974; George 
and Kim, 1976). The spectrum can be written as 


3 2 5 

$ (k) - 64 n v A 

ww' 


, 2 , 2 

k^ + k£ 


(1 + 4w 2 A 2 k 2 ) 3 


(3.12) 


where A is the turbulence length scale in m, and w is the turbulent 
intensity in m/sec. Although it is fairly accurate at low 
frequencies, the Dryden spectrum is not a very good model for high 


frequencies . 
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The Karman spectrum has been used by Amiet in his inflow 
turbulence analysis (1976). This model represents a more realistic 
energy distribution at high frequencies. With some rearrangements, 
the von Karman spectrum can be written in a form similar to the Dryden 
spectrum, i.e. equation (3.12): 


ww 


00 - 


55 

36 ^-^ 


£0/3) 

r(5/6) 


w 2 A 5 



r(i/3) 12 

r(5/6)/^rJ 


A 2 k 2 ) 


17/6 


(3.13) 


or, one can plug in the numerical values for Gamma functions and get 


* w (k) “ 0-^976 w 2 A 5 


(1 + 


k 2 k 2 

K 1 + *2 

1.7929A 2 k 2 ) 17 / 6 


Figure 3.6 shows a comparison of the two spectra. At high 
frequency, i.e. k » 1, the Dryden spectrum has a k'^ behavior, but 
the von Karman spectrum decays with a slower rate (k"^ 7 ^). At low 
frequencies, i.e. k « 1, although both spectra have the same k 2 
slope, their absolute values are quite different. 

As will be shown later in this chapter, for full scale helicopter 
rotors and wind turbines, the Karman spectrum yields better 
comparisons between the experiments and the predictions than the 
Dryden spectrum. On the other hand, for small scale tests, namely 
those tests involving model rotor testing in wind tunnels, 
calculations made with the Dryden spectrum are more accurate than 
those using the Karman spectrum. More attention has to be devoted to 
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turbulence spectral modeling equations to clear up the uncertainty of 
which model should be used. 

Wh«n turbulence is drawn to the rotor disk, severe distortion can 
occur to the ingesting flow. Amiet has studied the change of 
turbulent spectrum due to the mean flow contraction using a rapid 
distortion theory (Simonich et al . , 1986a, 1986b; Amiet et al., 

1986b). The mean flow contraction produces a nonisotropic turbulence 
field at the rotor disk. For hover and low speed vertical ascent, the 
inflow distortion results in an increased noise of 10 dB at low 
frequencies, where the acoustic spectrum is governed by the 
peak -and -valley structure, and an increase of 5 dB in high frequency 
broadband noise. However, at high speed forward flight, the mean flow 
distortion does not exist, thus the resulting noise will not change 
relative to the noise calculated using an inflow of homogeneous 
isotropic turbulence (Amiet, 1986b). 

3 - 2 - 2 Numerical Techniques Associated with Karman Spectrum 

In order to implement the Karman spectrum in the method of George 
and Kim, several modifications had to be made to the analysis. 

Following the same procedure as George and Kim's original analysis 
(1-976) , the far field sound can be found to be 
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<S 1 (*,f)> - 0.7859 
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where 
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M z nc 

0 



k 2 - (w Q /U) 2 + k 2 
w Q - | f-nf)| 

To evaluate equation (3.14), one has to handle an infinite summation 
of an integral function. As shown by George and Kim (1976) , the 
infinite summation can indeed be trimmed to a finite number of terms 
Thus our major task is to find an efficient way to evaluate the 
integral in equation (3.14). 

Rewrite the integral in equation (3.14) as 

“ e 2 (J 2 (V>e) + J?(*e)) 

p 0 1 

dx -Z — - JJT 

0 f 1 + 0 1 f 1 + 1 . 7929A 2 © 2 1 

1 i . m2 7 1 > 


where 
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X 3 k 

P 


2 

f - nQ 


U 

e 2 * x 2 + p 2 


The above equation can be approximated by: 


A 

J n dx + 

o 


a - m 2 , 2 


380.392 M 2 c 2 A 14/3 


J 


dx 

<K 2 + P 2 ) 7 /3 


where A is a constant above which the asymptotic expressions are 
accurate to within a given limit, and IT represents the original 
integrand. 

Since the finite integral from 0 to A can be handle easily, we 
will focus our attention on the second integral (from A to infinity). 
This integral can be done by first finding the integration from 0 to 
infinitity and then subtracting the integral of the same integrand 
from 0 to A. The integration from 0 to infinity can be found 
analytically, which is 


f dx 

! <x 2 ; pv/3 


15 " r < 2 / 3 > , o 7 p-11/3 

(2P) 11/3 ru/3) ra/3) 


Thus the original integral can be found by evaluating the 


following equation: 
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A 

J n dx + 
0 


a • m 2 > 2 


380.392 M 2 c 2 A 14 / 3 


( 0.7P' 11 / 3 


I <3 ' 15> 
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3.2.3 Effects of Length Scale and Spectrum on Hoise 

The difference between the two turbulent spectral equations is 
quite apparent as shown in Figure 3.6. Although the Doppler shifts 
make it more complicated, one can roughly identify a given frequency 
radiation with the inverse time required for a blade to pass through a 
turbulence component of length 1/k, where k is the wavenumber. Thus, 
for frequencies on the order of 10 kHz, at a tip speed of 100 m/sec, 
one is interested in turbulence component wavenumbers of order 
100 m . For a full scale helicopter, the typical turbulent integral 
scale is of order 100 m, then the peak of the spectrum is at 
wavenumber k of order 0.01 m This implies that the high 

frequencies come from wave numbers 10^* higher than the inverse 
integral scale. Refering to Figure 3.6, we see that the difference 
between the two atmospheric turbulence models can be of order of 10 dB 
at these wavenumbers . 

From another perspective, Figure 3.7 shows comparisons between 
von Karman and Dryden inflow turbulence noise calculations for a full 
scale UH-1 helicopter main rotor for both 0.1 m and 67.0 m integral 
scales. It is apparent that for small integral scales and low 
frequencies the inflow turbulence noise results for different 
turbulent models are in close agreement, but the difference becomes 
marked at high frequencies. 
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3.3 BOUNDARY LAYER/TRAILING EDGE NOISE 

Previous analyses of the boundary layer/trailing edge noise 
assumed zero blade angle of attack (Kim and George. 1980; Schlinker 
and Amiet. 1980). In practice, to produce desired loadings, rotor 
blades are seldom operated in such condition. Rather, they are 
operated at various non- zero angles of attack. In this section, we 
will extend the analysis to account for the important effect of angle 

of attack on rotor trailing edge noise. 

We start from the result of Kim and George (1980) , which is given 
in equation (3.9). Examining this equation, we can see that a change 
in angle of attack will affect the noise radiation through the changes 
on S the surface pressure spectrum. With the normalization factor 

pp- 

suggested by previous investigators, S pp can be written in the 
following form 

s pp ■ q2 r s o (S) °- 16) 

In order to accurately predict the trailing edge noise, clearly 
we have to find better approximations for both S , the displacement 
thickness, and the normalized pressure spectrum S Q (w) . 

3.3.1 Displacement Thickness 

The displacement thickness (5*) on rotor blades can be influenced 
by a number of parameters such as Reynolds number, Mach number, blade 
angle of attack, etc. Previous studies used flat plate boundary layer 
theory to calculate 6* and used the result as an input to the 
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analysis. As pointed out by Schlinker and Amiet (1980), the plate 
boundary layer theory cannot predict the correct displacement 
thickness for a rotor blade, except approximately for the zero lift 
case. Theoretically, displacement thickness is affected by Reynolds 
number, Mach number, angle of attack, etc. With increasing Reynolds 
number, 5* decreases slowly with Re* 1 / 5 ; with increasing Mach number, 
the compressibility effect tends to increase 5*. 

Brooks and Marcolini (1984a) made a number of measurements of the 
displacement thickness on two-dimensional airfoils of various size. 
Their results suggest that after appropriate normalization, the 
displacement thickness can be expressed in functional form as 

6* - 6o(Re,M) F(a) (3 1?) 

"A" 

where S Q is the displacement thickness on an airfoil with zero angle 
of attack, and F(a) is the correction factor for changes on angle of 
attack . 

The displacement thickness at zero angle of attack, after 
normalization by the airfoil chord, can be expressed as a function of 
Reynolds number, Re, and Mach number, M. In the present study, we use 
the experimental correlation given by Bies (1966): 



(3.18) 
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To find F(a) , we use the experimental results of Brooks and 
Marcolini (1984a). A curve fitting technique led to the following 
empirical expression for F(a) : 

F(a) - 1 + 0.028a, for a < 0° , and 

F( a ) - 1 + 0.028a + 0.016942a 2 ' 56 , for a > 0° 

For an airfoil with positive angle of attack, the first equation gives 
the correction factor of the displacement thickness at the lower 
surface and the second equation gives the correction factor at the 
upper surface. The above correlations are generally valid for a 
within -20° and 15°, and for Re within 8 x 10 4 and 5 x 10 6 , and for 
M < 0.25, where the boundary layer is turbulent and attached. This 
curve is shown in Figure 3.8. 


3.3.2 Surface Pressure Spectrum 

Next we examine S pp , the spectral density of the fluctuating 
surface pressures. As can be seen in equation (3.17), the term that 
is still left undetermined is S Q (w). An empirical expression for 
S 0 («) can be obtained from experiments. In the present study, two 
sets of experiments were used, by Yu and Joshi (1979), and by Brooks 
and Hodgson (1980). Their data agree fairly well with each other, and 
curve fitting with the basic high frequency form of the spectrum leads 
to the following expression (Chou and George, 1983). 


s 0 <“) “ 


1 . 732x10" 


(1 - 5.489w + 36 . 74w 2 + 0.1505w 5 ) 


(3.19) 


for 0.01 < u < 0.06, (w 


2xf6*/U) , and 
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a (DEGREE) 


Figure 3.8: Plot of F(a) vs. a. 

Sq ( w) - - - 1-4216x10 3 * U) 

(0.3621 + 4.1837w + 22.8185^ + 0.00135 3 + 0.00285 5 ) (3-20) 

for 0.06 < 5 < 20. 

Figure 3.9 shows the plot of Sq(co) vs. u> along with the range of 
experimental data and the flat plate result. 

3 - 3 - 3 Effects of Blade Angle of Attack on Broadband Noise 

Figure 3.10 shows the effect of angle of attack on rotor trailing 

edge noise demonstrated by calculations for a UH-1 helicopter main 
rotor. The results lead to the conclusion that the primary effect of 
an gl e °f attack is in the low to mid frequency range, where the noise 
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Figure 3 . 9 : Plots of Sq(co) vs. u> with Equations ( 3 . 19 ) and ( 3 . 20 ). 

level increases with angle of attack. However, in high frequency 
range, this change of noise level due to change of angle of attack is 
not very significant. The comparison of predictions using present 
analysis to that using flat plate data only (e.g. Kim and George, 

1980) shows the important effect of the angle of attack on rotor 
trailing edge noise. 

Trailing edge noise is only one source of rotor broadband noise. 
Other mechanisms such as inflow turbulence and tip vortex separation 
also contribute significantly to the noise radiation. Trailing edge 
noise can be important for low inflow turbulence cases or when 
considering a large sized rotor (George and Chou, 1983a, 1984a). Thus 
to evaluate the present analysis by comparing with existing 

t also include other possible sources. 


experiments, one mus 




FREQUENCY, Hz 


Figure 3.10: Effect of Angle of Attack on Rotor Trailing Edge Noise 

UH-1 Main Rotor, 4 , - -27°. 
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3.4 TIP VORTEX FORMATION NOISE 

The original analysis of turbulence/trailing edge interaction in 
the tip vortex separation area used a simplified cross-flow analogy to 
estimate the turbulence level and the surface pressure spectrum in the 
tip separation region (George et al . , 1980). In the present study, 
this analysis is extended, with the pressure spectrum being replaced 
by that associated with the leading edge vortices on delta wings 
(George and Chou, 1984a). 

3.4.1 Surface Pressure Spectrum from Delta Wings 

Again we start with the result of George et al . (1980), i.e. 
equation (3.11). As mentioned previously, the accuracy of the 
prediction depends upon the accuracy of the model used for describing 
the surface pressure fluctuation and upon parameters such as the 
extent of the tip separation region. 

The original analysis used surface pressures measured from the 
separation region behind a 2-dimensional fence. The model used in the 
present study is based on the pressure spectra measured under the 
similar edge separation vortices on delta wings. As sketched in 
Figure 3.2, the delta wing flow geometry is very similar to helicopter 
rotor tip flow, including both the primary and secondary separations 
from the edge and the axial outer flow. The separation geometry is 
also influenced by rounded or sharp edges as in the delta wing tip 
case. Our goal is to construct a suitable correlation for the 
pressure spectra in the delta wing case and relate this correlation to 
the rotor tip flow case of interest. 
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Richard and Fahy (1965) have analyzed the turbulent flow beneath 
the leading edge vortices of several delta wings of different 
planfornrs and various angles of attack. They presented spectra from a 
number of investigators, non-dimensionalized in several ways, however, 
none of which were satisfactory for our case. 

In order to find a normalization suitable for application to our 
tip vortex separation case, we first studied the delta wing flow 
geometry and pressure data measured in the comprehensive experiments 
of Peckham (1961). Based on flow visualization results, the locations 
of peak negative pressures, and on pressure distributions, the value 
of the transverse separation scale L was found to be related to the 
local chord length. The edge shape is definitely important as noted 
also by Barlett and Vidal (1955). 

Next the maximum negative pressure coefficients under the 
v o^ti.ces relative to those on the nearby surface were used with the 
Bernoulli's equation to estimate the maximum velocity ratio V /U. 
Assuming that the velocity in the vicinity of the vortex is 
approximately the same as that on nearby surface, one obtains: 


v m /u 

nr 


^- c p . > 

r min 


1/2 


Then we used the estimated values of V m and L to normalize the 
spectra given by Richard and Fahy (1965). The results, as shown in 
Figure 3.11, give a fairly good correlation. 

With the usual curve fitting techniques, these spectra can be 
approximated by 

log 10 (S 1 (w)) - -3.475 - 1.654(o + 0.82) 2 


(3.21a) 
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(3.21b) 


for a < -0.82 (o ■ log-^Qw) , and 

lo gl 0 ( s i(w)) - -3.475 - 0.984 (a + 0.82) 2 

for a > -0.82. But S 1 (5) - 0 for 5 < 0. The fitted curve is plotted 
against data as shown in Figure 3.11. 


3.4.2 Other Parameters 

In order to define the appropriate spectrum for a rotor tip case, 
we need the estimates of L and V m for rotor tip flows at different 
angles of attack. The flow visualization, pressure measurements, and 
velocity measurements of Gray et al . (1980) and of Chigier and 

Corsiglia (1971) were used to estimate values of L/c and V m /U. 

Figures 3.12 and 3.13 show the plots of L/c and V m /U versus a 
respectively. V m /U can be expressed as 


V /U * 1 + 0.0359 a 
nr 

where a Is the tip angle of attack in degrees. 

The spanwise extent of the separation bubble remains to be 
determined. George and Chou (1984a) suggested the following 
correlations 


(3.22) 


L/c - 0.023 + 0.0089 a 


for square tip blades, or 


(3.23a) 


L/c - 0.0074 (a - 2) , for a > 2 C 


(3.23b) 


for blades with rounded tips. 
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Figure 3.12: Maximum Circumferential Velocity Ratio vs. Angle of 

Attack . 

Brooks and Marcolini (1984b) have performed a series of 
experiments on the tip vortex formation noise from airfoils with 
rounded tips. They found that the basic scaling of the tip vortex 
formation noise appears to be correct. The correlation for V m as 
given in equation (3.22) is found to be consistent with the 
experiments (within 5 to 10 percent). However, they also found that 
if one chooses 4 to 5 percent turbulent intensity contours as a guide 
to determine the separation length scale L, the correlation given in 
equation (3.23b) does not seem to be appropriate. They suggested that 
equation (3.23b) should be replaced with 


L/c — 0 . 008 a 


(3.24) 
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Figure 3.13: Lateral Extent of the Separation Bubble on Rotor Tips 

vs. Angle of Attack. 

3 . 4.3 

It is clear that tip vortex formation noise is extremely 
sensitive to the change of angle of attack and the shape of the rotor 
tip. Figure 3.14 shows the results of some calculations for a UH-1 
helicopter main rotor. Calculated spectra using (3.23) for L are 
shown for pitch angles of 10° and 15° for both square and rounded tip 
rotor. Also shown in the figure are the corresponding boundary 
layer/trailing edge noise calculations. Clearly, tip vortex formation 

e-a 
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noise is favored by high angles of attack, square tips, and wide tip 
chords (low aspect ratio, untapered blades). 

3.5 TRAILING EDGE THICKNESS NOISE 

Turbulent vortex shedding from blunt trailing edges is a source 
of rotor high frequency broadband noise (Schlinker and Brooks, 1982) . 
Brooks and Hodgson (1980) were the first to demonstrate the importance 
of this self - generated noise mechanism. In their experiment with an 
isolated airfoil, they found that the radiated noise increased 
significantly for airfoils with thickened trailing edges. 

The first attempt to study this noise mechanism theoretically for 
rotors was done by Chou and George (1985). They studied the 
parametric dependence of this noise mechanism and developed a method 
to predict rotor broadband noise associated with turbulent vortex 
shedding from blunt trailing edges. Further work is needed to 
definitively establish the turbulent vortex-shedding properties of 
wings and rotor tips. 

3.5.1 Surface Pressure Spectrun 

Blunt trailing edge noise radiation is a result of higher surface 
pressure fluctuations near the airfoil's trailing edge due to 
turbulent vortex shedding. To predict this noise mechanism, our first 
task is to scale the fluctuating surface pressures. Using dimensional 
analysis or physical reasoning, the following relationships for the 
parametric dependence for Spp and f were found 

S - Q 2 t 3 / U 

PP 


f - U / t 
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Figure 3.14: Effect of Angle of Attack on Tip Vortex Formation 

Noise, <f> - -27°, UH-1 Main Rotor. 
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where S is the power spectral density of surface pressure 
PP 

fluctuation, f is the frequency in Hz , t is the trailing edge 
thickness and U is the free stream velocity. By using the above 
relationships, the surface pressure spectrum S pp is assumed to have 


the following form 


S pp (f) - q 2 t 3 S 2 (S) / U 


(3.25) 


where u * 2nft/ U. The normalized spectrum S 2 (w) can be found 
empirically from experimental data. Due to the lack of measurements 
for surface pressure fluctuation near blunt trailing edges, such data 
were obtained by inverting the acoustic data measured by Brooks and 
Hodgson (1980) using their stationary airfoil analysis. Using the 
above scaling relation, the data representing a wide range of free 
stream velocities and trailing edge thickness are collapsed reasonably 
well to a single curve. By using curve fitting techniques, the 
empirical expression for the normalized spectrum S 2 is found as 
follows 


log 10 S 2 (S) - 17.394 - 106.57a - 158.12a 2 + 99.27a 3 
-33.249a 4 + 16.721a 5 

for 0 . 2 < a < 2 . (a- log 1Q (w)), otherwise 


S-(«) “ 0 


(3.26a) 


(3.26b) 


Figure 3.15 shows the fitted curve along with the experimental data. 

To predict the additional broadband noise radiation from a rotor 
due to its blunt trailing edges, several assumptions have to be made. 
First, the source is modelled as rotating radiating dipoles; this 
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1 0 * LOG { SppU/q^t 3 ) 



Figure 3.15: Normalized Surface Pressure Spectrum S 2 («) and the 

Experimental data of Brooks and Hodgson (1980). 

assumption can be justified from the acoustic data measured by Brooks 
and Hodgson (1980), in which a strong dipole radiation pattern was 
observed. Then the noise can be predicted with an analysis similar to 
that of Kim and George for boundary layer/trailing edge noise (1980), 
with surface pressure spectrum replaced with (3.26). The result of 
this analysis is given in equation (3.9). 
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3.5.2 Effect of Trailing Edge Thickness on Broadband Noise 

Figure 3.16 shows the effect of trailing edge thickness on rotor 
broadband noise. Calculations were made based on the low speed fan as 
used in Lowson's experiments (1972). It is clear that trailing edge 
thickness is a very important parameter for the rotor noise problem. 
Generally speaking, the noise spectra due to turbulent vortex shedding 
from blunt trailing edges are peaks occurring at various frequency 
ranges. The peak frequency depends on the trailing edge thickness, a 
small trailing edge thickness will generate a high frequency peak and 
a thick trailing edge will result in a peak of lower frequencies. The 
level of the spectrum peak also depends on the thickness of trailing 
edge; peak level increases roughly according to the third power of 
trailing edge thickness. 

In conclusion, turbulent vortex- shedding noise from blunt 
trailing edges is a very important broadband noise source for rotors. 

A slightly blunted rotor trailing edge can contribute significantly to 
the overall noise spectrum. The present analysis provides reasonable 
predictions for such mechanisms. However, a more accurate prediction 
could be achieved with a better empirical expression for S 2 (5) , the 
normalized spectrum. Clearly, more measurements of surface pressure 
fluctuations near blunt trailing edges are needed. 

3.6 RESULTS AND DISCUSSIONS 

In the present work, we reviewed and made significant extensions 
to rotor broadband noise analyses to allow more accurate prediction of 
rotor noise spectra. Our analyses, although evaluated by computer 
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programs, are primarily analytical and thus helpful in understanding 
the nature of the noise generation. The sources considered are load 
or surface pressure fluctuations due to inflow turbulence, due to 
turbulent boundary layers passing the blades' trailing edges, due to 
tip vortex formation, and due to turbulent vortex shedding from blunt 
trailing edges. 

Comparisons to more numerically based approaches show that our 
analyses are accurate but restricted to advance ratios less than 
approximately 0.4 (which include all cases of practical interest). 

The present study neglects in-plane noise mechanisms and is thus 
restricted to angles which are not very close to the rotor plane. 

The results of the present study indicate that inflow turbulence 
noise depends strongly upon ambient conditions and dominates at low 
frequencies. Trailing edge noise and tip vortex noise are important 
at higher frequencies if inflow turbulence is weak. Boundary layer 
trailing edge noise is important, especially for large sized rotors. 
This noise increases slowly with angle of attack but not as rapidly as 
tip vortex noise, which can be important at high angles of attack for 
wide chord, square tips rotors. Turbulent vortex shedding noise from 
blunt trailing edges can be very important in mid- to high-frequencies 
depending upon the type of rotor and its trailing edge thickness. 

3.6.1 Comparisons of Analysis to Experiments 

An extensive search was made of existing experiments and 
calculations based on the various prediction methods were made for 
comparison purposes. This study shows that present analyses are 
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adequate to predict the spectra from a wide variety of experiments on 
fans, full scale and model scale helicopter rotors, wind turbines to 
w ithin. about 5 to 10 dB. Better knowledge of the inflow turbulence 
improves the accuracy of the predictions. 

^ comparisons, the data estimated as input to the analyses 
and the correlations are given in the figure captions. Other input 

were taken from the experimental papers. Inflow velocities 
were estimated using simple momentum theory with thrusts determined by 
simple blade element theory. As mentioned above, the inflow 
turbulence values were often estimated. In cases where separate 
calculations are shown for separate mechanisms, the results should be 
summed in order to compare to the experiments. However, in order not 
to clutter the figures, this was not done in most cases. 

First we compare the present analyses to the experimental data 
for full scale helicopter rotors. Generally there are two types of 
full-scale helicopter rotor tests available: the whirl tower test and 
the flight test. The whirl tower test has several advantages over the 
flight test. Since only one rotor is involved, there is no problem 
with aerodynamic interactions with other rotors such as main 
rotor/tail rotor interactions. Also, other polluting noises such as 
noises generated by drive motor and gear box, etc. are comparatively 
easy to control. Therefore, these tests are considered cleaner than 
the flight tests. However, the flight tests do give more information 
on the overall helicopter noise problem. 

First we compare the present analyses to a whirl tower test by 
Leverton (1973). He tested a full scale S-55 rotor on a test rig in 
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an inverted position in order to eliminate the effect of recirculation 
which occurs when a rotor wake is directed toward the ground. The 
spectra measured were taken from a tethered balloon at various angles 
to the test rotor plane. His tests varied both load and RPM. The 
primary missing information in Leverton's results is any data on the 
inflow turbulence. Neither the turbulent intensity nor the scale were 
measured. As the inflow was drawn from near the ground, the turbulent 
integral scale could be quite reasonably estimated from the fairly 
well established empirical relationship that 

A - 0.9 h 

where A is the turbulent integral scale and h is the height above the 
ground (Etkin, 1961). Similarly, values for turbulent intensity for 
various weather conditions can also be estimated from the extensive 
data and correlations in Lumley and Panofsky's monograph (1964). 

The first set of data we choose were taken at an angle of -75 
from the rotor plane where all of the analyses would be expected to be 
within the range of their assumptions. Figure 3.17 shows the 
comparison of a range of predictions. It is also clear that at the 
lower frequencies, say below 1000 Hz, the boundary layer trailing edge 
noise and the tip vortex noise mechanisms are not important. However, 
at frequencies above 1000 Hz they become quite important, with 
boundary layer trailing edge noise being the more important one in 
this case. Fink's boundary layer noise correlation is seen to be a 
reasonable approximation to the more exact boundary layer trailing 
edge noise calculations. Most of the noise below 1000 Hz is shown to 



FREQUENCY, Hz 


Figure 3.17: Comparison of the Range of Predictions for a Full Scale 

Helicopter Rotor with the Experimental data of Leverton 
(1973), 4 > - -75°, A - 0.57 m, w - 1 m/sec. 




-97- 


be inflow turbulence noise based on the estimated turbulent 
properties. Both the analyses of George and Kim (1976) and of Amiet 
(1976) agree within 5 dB with each other and with Leverton's data. 

In Figure 3.18, a comparison is shown among calculations based on 
the three mechanisms of George and co-workers, the two of Amiet, and 
some data of Leverton at an angle of -11.5° from the rotor plane. As 
all of these analyses ignore in-plane force components and as George 
and co-workers use a blade dipole directivity, the agreement would not 
be expected to be quite as good. It is not clear, however, how many 
of the differences are due to which of these effects. The inflow 
turbulence no 

however seem to be in better agreement with the experiments in this 
case . 

Next we examine the analyses against the experiments of full 
scale wind turbines. Figure 3.19 shows the comparison to full scale 
MOD-2 wind turbine data presented by Hubbard et al. (1981). The 
background noise was measured and shown to be well below the measured 
spectra. But again no information on turbulent intensity or integral 
length scale was given. Thus the calculations are based on estimated 
turbulence properties, aided slightly by the fact that at least the 
wind speed is known. There is little question that the predictions 
using the Karman spectrum are in better agreement with the experiments 
than those using the Dryden spectrum. The results shown used 
calculated trailing edge thickness, which is based on standard NACA 
23XXX series airfoil sections as in the MOD-2 blade application. The 
results compare favorably to the acoustic data obtained from the 



FREQUENCY, Hz 


Figure 3.18: Comparison of Predictions with the Experimental Data of 

Leverton (1973), 4 > - -11.5°, Other Parameters are Same 
as in Figure 3.17. 
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experiment. It is clear that the primary noise sources above a few 
hundred Hertz are boundary layer/trailing edge noise and trailing edge 
thickness noise. 

The noise from small low speed axial flow fans is also 
considered. The indoor experiments of Lowson et al . (1972) are chosen 

for our comparison. Tests were run in an anechoic room for both a 
ducted and an unducted fan, both before and after recirculation built 
up in the room. The non- recirculation cases are chosen for our 
comparisons since test conditions are better defined. RPM, tip angle, 
and tip shape were varied as well. The turbulence was measured by a 
hot wire anenometer with limited frequency response in this 
experiment. We estimated the turbulent integral scale as 0.1 meter. 
Figure 3.20 shows the comparison of the present analyses with Lowson' s 
experimental data. Again all possible broadband noise sources are 
included. The results show excellent agreement with the experiment. 

The model scale helicopter rotor tests are also considered. A 
set of high quality tests has been carried out by Paterson and Amiet 
in the UTRC anechoic wind tunnel facility (1979). In these tests, 
both vertical ascent and forward flight were simulated and different 
grids were used upstream to generate controlled inflow turbulence. 
Measurements were made of both the turbulent intensity and integral 
scale so that, in these cases, none of the parameters needed to be 
estimated. Figures 3.21 through 3.23 show comparisons of calculations 
to data presented by Paterson and Amiet. In the no grid case (low 
inflow turbulence), it is clear that both tip vortex and boundary 
layer noise are important at the higher frequencies. In all the 
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cases, the calculations and experiments agree to within about 5 dB . 
Unlike the full scale rotor cases, inflow turbulence noise 
calculations using the Dryden spectrum agree better with the 
experiments than those using the Karman spectrum. In Paterson and 
Amiet s original report, they had also shown good agreement with 
Amiet's inflow turbulence noise analysis whenever the primary noise 
source was inflow turbulence. 

3.6.2 Comparisons of Analysis to Each Other 

In this section the results calculated by the methods of Amiet 
and of George and co-workers are compared to show the effects of 
different assumptions in the analyses. The computational approach of 
Amiet allows treatment of forward flight (non-zero advance ratios) and 
more accurate basic blade noise directivity. The George and Kim 
approach has been implemented for both the Von Karman model and the 
Dryden model of the inflow turbulence spectrum. We will examine each 
of these effects by comparing the results of the calculations made by 
different methods. 

Figure 3.24 shows the effect of forward flight on inflow 
turbulence noise as calculated by Amiet's analysis and compared to 
hover calculations based on George and Kim's analysis. It is notable 
that the advance ratio effect is not very important for any case of 
interest to helicopters (i.e. advance ratio below 0.4). 

Similarly, Figure 3.25 compares changing advance ratio for 
boundary layer/trailing edge noise. Here, the basic inputs vary since 
the calculation using the Kim and George approach uses an airfoil 
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FREQUENCY, Hz 


Figure 3.24: Effect of Advance Ratio on Rotor Inflow Turbulence Noise 

Calculation, UH-1 Main Rotor, 4 > - -60°, A - 67 m, 
w - 1 m/sec . 
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boundary layer thickness correlation rather than the flat plate 
results incorporated in their publications and in those of Amiet. In 
this boundary layer trailing edge noise case, the results again show 
that the effects of advance ratio are not important for values less 
than 0.4. 

As discussed in a previous section, Amiet' s computational model 
incorporates an accurate basic blade noise radiation directivity for 
the pressure normal to the blade mean chord line. The methods of 
George and co-workers approximate the basic directivity by a dipole 
normal to the rotor plane, which would be expected to lead to 
underestimates for angles near the rotor plane. Both Amiet' s and 
George and co-workers' analyses ignore in-plane forces/noise 
mechanisms. Figures 3.26 through 3.28 compare the directivities for 
both inflow turbulence and trailing edge noises. 
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Figure 3.25: Effect of Advance Ratio on Rotor Trailing Edge Noise 

Calculation, UH-1 Main Rotor, <f> - -90°. 
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Figure 3.26: Directivity of Rotor Inflow Turbulence Noise 1500 Hz and 

6000 Hz. S - 55 Full Scale Rotor, A - 0.57 m, w - 1 m/sec. 
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Figure 3.27: Directivity of Rotor Trailing Edge Noise, 1500 Hz, S- 

Full Scale Rotor. 
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Figure 3.28: Directivity of Rotor Trailing Edge Noise, 6000 Hz, S-55 

Full Scale Rotor. 



Chapter IV 

HELICOPTER TAIL ROTOR NOISE 


In this chapter, various helicopter tail rotor harmonic arid 
broadband noise mechanisms are studied and new source models are 
proposed. Of particular importance are those sources due to 
interactions with the wakes of the main rotor, the main rotor hub, the 
fuselage,, and the engine exhaust are examined. The flow field around 
the tail rotor and the noise it generates are modeled using available 
aerodynamic and acoustic analysis along with some necessary estimates 
of turbulence properties. Representative calculations show that main 
rotor wake is the strongest contributor to both tail rotor harmonic 
and broadband noise. The fuselage separation wake, the engine 
exhaust, the main rotor hub, and the hub -shaft wakes are important to 
the tail rotor broadband noise only. The main rotor tip vortices are 
not important to tail rotor broadband noise. However, they have a 
major effect on the tail rotor harmonic noise; details of the tail 
rotor blade-main rotor tip vortices interaction will be discussed 
later in chapter V. The goal of the study discussed in this chapter 
is to determine the relative importance of various tail rotor harmonic 
and broadband noise mechanisms. 
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4.1 REVIEW OF PREVIOUS TAIL ROTOR NOISE RESEARCH 

As mentioned in the Introduction, the importance of tail rotor 
noise was well appreciated by the year 1970 (Lynn, 1970). The 
complexity of the flow field around tail rotor has also been studied 
for some time because of its impact on tail rotor sizing, blade stall, 
fin interference, and directional control of helicopters. Some of the 
flows which have been identified are: 

1) interaction of the tail rotor with the mean downwash flow field 
of the main rotor; 

2) interaction with the rolled-up "ground vortices"; 

3) interactions with the main rotor blades' individual tip 
vortices; 

4) interactions with the main rotor disk's "wingtip vortices" (from 
the overall horseshoe vortex due to the mean lift of the main 
rotor disk) ; 

5) interactions with the wakes of stabilizer, fin, and tail boom; 

6) operation of tail rotor in the "vortex ring" state. 

One important finding of these studies was that it is desirable 
for the tail rotor to rotate "bottom forward" to minimize interactions 
with ground and wingtip vortices. Later, it was realized through 
other experiments that this direction of rotation is also desirable 
for noise reduction (Leverton, 1977, 1980). This is because the noise 
due to high speed interactions of the tail rotor blades with the main 
rotor tip vortices seem to be beamed in the blade motion direction. 
Leverton (1980) reported noise reductions of approximately 15 dBA in 
approach and 5 dBA overhead for an optimized tail rotor design. 
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While the transonic, near parallel and streamwise interaction 
characteristic of main rotors has been treated by several 
investigators (George and Chang, 1983c, 1984b; George and Lyrintzis, 
1986; Srinivasan and McCroskey, 1984; McCroskey, Yu, and Smetana, 
1984), the high speed, near normal incidence blade-vortex interaction 
characteristic of tail rotors has not been studied to date. The low 
speed, near normal incidence case has been studied by Schlinker and 
Amiet (1983), Amiet (1986), and by Ahmadi (1984a, 1984b). 

The importance of wake ingestion on rotor- generated noise was 
shown in the study of Levine (1976). He reported increases of 5 to 10 
dB in both narrowband and broadband noise when a Sikorsky S-58T 
operated with the main rotor wake being blown into the tail rotor. 

Another experimental study of tail rotor-main rotor wake 
interaction noise involved wind tunnel tests of a model scale UH-1 
with tail rotor position and direction of rotation as test parameters 
(Balcerak, 1976; Pegg and Shidler, 1978). Balcerak (1976) made a 
parametric study which varied tail rotor location, fin blockage area, 
tail rotor rotation direction, rotor speed, thrust, and the tail rotor 
operating mode (i.e, pusher/tractor). Later, Pegg and Shidler (1978) 
tested the same model, extending the work and emphasizing 
identification of the aeroacoustic mechanisms producing the noise. 

They found an approximately 12 dB increase in broadband noise when the 
main rotor flow was drawn through the tail rotor and significant 
increases in harmonics under various conditions. Combination tones 
due to the main rotor shed vortices and turbulent wake were apparent 


in the non-zero advance ratio cases. 
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Some of the representative results of Balcerak (1976) are given 
in Figures 4.1 through 4.4. Figures 4.1 and 4.2 give the noise 
spectra for isolated main and tail rotors. Figure 4.3 shows the 
measured noise spectrum in a hover condition when both rotors are 
operated. Notice that it is essentially the sum of the spectra shown 
in Figures 4.1 and 4.2. Figure 4.4 shows the noise spectrum measured 
for a forward flight condition (advance ratio based on main rotor tip 
speed is 0.28) when both rotors are operated; significant increases in 
both harmonic and broadband noise are apparent. These studies also 
showed, as did Leverton's (1977, 1980), that interaction noise is more 
severe when the advancing rotor blades interact with the main rotor 
wake. Further discussions of these and Leverton's experiments were 
given by White (1978) . 



Figure 4.1: Noise Spectrum, Main Rotor Only, from Balcerak (1976). 
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Figure 4.2: Noise Spectrum, Tail Rotor Only, from Balcerak (1976) 
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Figure 4.3: Noise Spectrum, Both Rotors, Hover, from Balcerak (1976) 
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Figure 4.4: Noise Spectrum, Both Rotors, p - 0.28, from Balcerak 

(1976). 


Tail rotor noise reduction techniques have been studied by Levine 
(1976) using the methodology available at that time. He modeled the 
wake using a simple average rigid wake model of Heyson (1961). This 
model assumes the wake to be a skewed cylinder of radially varying 
vorticity of the same diameter as the main rotor. Using the 
velocities from Heyson' s model, Levine computed the low-order airload 
harmonics on the tail rotor blades. Then, the noise due to the 
azimuthally- varying blade loading is calculated using the analysis of 
Lowson and Ollerhead (1969). This approach obviously leaves out most 
of the complexities of a real problem. 
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4.2 PRESENT RESEARCH 

In the present study, sample calculations are made for a UH-1 
type helicopter in order to compare the noise spectra from the various 
mechanisms. In each case the helicopter is assumed to fly in a 
condition such that maximum interaction occurs. For example, when 
calculating fuselage wake effects, the full wake width is assumed to 
impinge on the tail rotor. 

4.2.1 Isolated Tail Rotor Noise 

In the present study, isolated tail rotor broadband noise is 
calculated using the analyses we discussed in chapter III. The 
results are compared to those from the main rotor to determine the 
relative importance of isolated tail rotor broadband noise. 

Mechanisms considered are inflow turbulence noise (due to the 
interaction with ambient turbulence), boundary layer/trailing edge 
noise, tip vortex formation noise, and trailing edge thickness noise 
(due to turbulent vortex shedding) 

4-2.2 Harmonic Noise Due to Non -uni fora Inflow 

In order to predict the tail rotor harmonic noise due to the 
interactions with various non-turbulent mean wakes, three steps are 
necessary : 

1) Predict the inflow velocities at the tail rotor disk (which 

include the wakes of main rotor, main rotor hub/shaft, fuselage, 
and engine exhaust) ; 
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2) Calculate tail rotor blade airload harmonics using the known 
incident velocity; 

3) Calculate the harmonic noise from the known loading harmonics. 

The comprehensive rotorcraft prediction code (CAMRAD) of Johnson 
(1980, 1981) is used for the prediction of main rotor mean wake. This 
code calculates the main rotor blade section loadings based on 
non-uniform induced inflow and on the forward- flight free-wake model 
of Scully (1975). The induced velocities at the tail rotor disk can 
then be found. Figure 4.5 shows the sketches of the boundaries of the 
main rotor wake at various air speeds. 


90 knot 5 



Figure 4.5: Sketches of the Main Rotor Wake and Engine Exhaust Flow, 

from Leverton (1977). 

For tail rotors with rotational speeds which are different from 
that of the main rotors, the CAMRAD code does not allow the detailed 
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loading calculation on individual tail rotor blades. Instead the code 
uses the velocity at the tail rotor hub to find the tail rotor blades' 
loading. Thus, the noise calculated from these loadings may be 
inaccurate if the interactions with the smaller scale velocity 
perturbations in the main rotor wake are dominant. 

The engine exhaust mean flow is approximated by the far field 
results of turbulent jets (Schlichting, 1979): 


3K 

87re 0 x(l+r7 2 /4) 2 

(3K) 1/2 (qV/4) 
r “ 4* 1 /2 x(1+ „2 /4) 2 


where 

K 

b 


1/2 


- ( 159b l/2 U jet> 

- radius of the engine exhaust jet (- 0.848x) 

- ((3K) 1/2 /(4 ? t 1 / 2 € q ) ) (y/x) 

- 0.0256b 1/2 lh ec 

- streamwise coordinate 


radial coordinate 


(4.1) 


For the fuselage wake, we estimated the velocity field at the 
tail rotor using a velocity distribution: 



C e 


- (y 2 /^ 2 ) 


(4.2) 
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The values of C and ft are estimated from the measurements on a BK-117 
helicopter (Huber and Polz , 1982) and scaled to UH-1D. The fuselage 
separation wake can also be similarly modelled by scaling the results 
of Polz and Quentin (1981) . 

For the main rotor hub, the characteristics of axisymmetric far 
wakes are: 



r c _^_ 1 1/3 - (y 2 /b 2 ) 

1 b 2 x 2 J 


where 

C d - drag coefficient 
b - (0.18C d Ax) 1/,:3 

A - frontal area of the rotor hub 


(4.3) 



Figure 4.6: Sketches of the Wakes from Fuselage and Main Rotor Hub, 

from Polz and Quentin (1981). 
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The loadings on the tail rotor are calculated with Johnson's 
CAMRAD code (Johnson, 1980, 1981) or with the quasi - steady , 
blade-element theory for rotors (Dommasch et al . , 1967). These 
loadings are then decomposed to find the harmonics using a Fast 
Fourier Transform algorithm (Bergland and Dolan, 1979). The harmonic 
noise is then calculated using the analysis of Lowson and Ollerhead 
(1969). This approach should be accurate enough for the present 
comparisons, but it would be better to use the non-compact methods of 
Farassat for more accurate results (Farassat, 1975, 1983). 

4.2.3 Broadband Noise Due to Various Turbulent Wake Effects 

The primary broadband noise mechanism is the inflow turbulence 
noise. Both Dryden and Karman spectrum are used in the present study 
to define the turbulent inflow. We need to estimate the turbulent 
intensities and the turbulent integral scales for the cases 
considered . 

The turbulence due to the overall main rotor wake was estimated 
in two ways : 

1) based on shear layer results for the edge of the wake, we 
estimate the intensity and the turbulent length scale as 
(Townsend, 1976) 

A - 0.5 (main rotor radius), and (4.4) 

w - 0.15 U ( ^ ownwas |_ i , or 

2) based on the blade drag, we estimate 

A - 0.4 (main rotor chord), and (4.5) 

w - 0.5 (main rotor blade speed). 



-123- 


These estimates are quite approximate, but there is currently no 
information available to support more accurate estimates. Some work 
has begun on turbulence in wind turbine wakes (Bossanyi, 1983) but it 
has not developed to the stage of being applicable to the rotor near 
wake region. 

In order to estimate the turbulence in the core of the main rotor 
tip vortices, several experimental studies are considered (Schlinker 
and Amiet, 1983; Biggers and Orloff, 1975; Tangier, 1978; and Owen, 
1970). From their results, it is estimated that 

A = vortex core diameter “ 0.15 (main rotor chord) , and 
w - °’ 13 U c , max 

where U is the maximum circumferential velocity of the main rotor 

c , max 

tip vortex. 

The turbulence properties of the engine exhaust jet and hub wakes 
were estimated based on the properties of similar flows ( Schlichting , 
1979; Tennekes and Lumley, 1972) . For the 2-D hub shaft wake, we 
estimate 


A =* hub shaft diameter, and 
w - 0.35 U 

where U is the forward flight speed of the helicopter. For the engine 
exhaust jet, we estimate 


A ** diameter of 


W - 0.25 U Jet 


where U. is the 
j e l 


the exhaust jet, and 


engine exhaust centerline 


velocity . 
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For the fuselage near wake, it is estimated that 

A =* 0.66 D , and 
sep 

w - 0.16 U 

where D sep is the diameter of the fuselage wake. 

The properties of the axially symmetric far wake are used to 
model the wakes of main rotor hub or the fuselage, we estimate 

A ~ D, and 
w - 0.2 U 


For an unseparated fuselage boundary layer we estimate 

A =* anc * 

w - 0.13 U 

where the boundary layer thickness on the fuselage. 

For each turbulent flow it is also necessary to estimate the 
fraction of the tail rotor disk immersed in the turbulence. For the 
present exploratory calculations, we assumed the maximum values 
possible for each case, i.e., the helicopter is trimmed so that the 
maximum amount of the flow in question (e.g., wake or jet) are 
incident on the tail rotor. 

For the engine exhaust jet, the effective tail rotor area is 
estimated to be 1/2 of the entire tail rotor disk; for the fuselage 
wake, 2/5; for the unseparated fuselage boundary layer, 1/8; and for 
the main rotor wake, the entire tail rotor disk is effective. For the 
tip vortices the intersection area was time-averaged over the passage 
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of the vortices through the disk; for ’'unburst’ 1 vortices the effective 
area is about 6 x 10"^ of the tail rotor disk, and for ’’burst” 
vortices the effective area is about 10 of the tail rotor disk. 

It should be noticed that these turbulent properties are not 
conservative estimates thus the results represent maximum noise that 
may be generated by the tail rotor. 

4.3 RESULTS AND DISCUSSIONS 

In the present study, calculations were made for parameters 
representative of a UH-1 helicopter. In order to calculate the 
interactions with the main rotor wake, conditions of level flight at 
approximately 100 mph are assumed. 

The first cases calculated were the basic broadband noise sources 
of the isolated main rotor and tail rotor, the results are shown in 
Figure 4.7. The sources include atmospheric inflow turbulence, 
boundary layer- trailing edge interaction, tip vortex 
turbulence - trail ing edge interaction, and trailing edge thickness 
noise. The last mechanism is calculated for both 1.27 mm (0.05 inches 
for typical UH-1 tail rotor blades) and for 2 . 54 mm trailing edge 
thickness. The results are given in terms of one Hertz bandwidth 
spectra for an observer at a distance of 62.5 meters to the rotor hub 
and at an angle of 63° from the tail rotor plane (-27° from the main 
rotor plane). It is clear that the isolated tail rotor generates far 
less broadband noise than the main rotor, except perhaps when the tail 
rotor blades' trailing edge are blunted. 
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Figure 4.7: Main Rotor and Tail Rotor Broadband Noise (Inflow 

Turbulence Noise Due to Ambient Turbulence) . 


Figure 4.8 presents the calculated tail rotor noise due to the 
interactions with the main rotor wake, which is calculated using 
Johnson's CAMRAD code. Clearly the lower harmonics are very 
important. The tail rotor broadband noise due to the main rotor wake 
turbulence is also very important; it can be on the order of 15 dB 
greater than the main rotor broadband noise in some frequency ranges, 
depending upon the actual turbulence scales. 

Figure 4.9 presents calculations for the tail rotor broadband 
noise due to interactions with the turbulence in the vortex cores of 
the main rotor wake. Clearly the noise is much less than that due to 
the main rotor turbulent wake. This is because the effective area in 
this case is much less than other tail rotor broadband noise 


mechanisms . 
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Figure 4.8: Tail Rotor Harmonic and Broadband Noise Due to 

b Interactions with the Main Rotor Wake. 
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Figure 4.9: Tail Rotor Harmonic/Broadband Noise Due to Interactions 

with Main Rotor Tip Vortices. 
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Figure 4.10 shows the tail rotor noise due to interactions with 
the fuselage wake. The separated flow harmonic noise is well below 
that calculated for the tail rotor - main rotor wake interaction shown 
in Figure 4.8. The broadband noise is quite important since the 
turbulent intensity is high and the effective area is large. The 
calculations using the analysis of Homicz and George (1974) show the 
peak and valley structure of the broadband noise at low frequencies. 
The attached fuselage boundary layer wake effects are also shown. 

Here the effect of turbulent length scale is evident; an attached 
boundary layer generates smaller scale turbulence in its wake compared 
to the fuselage separated flow case. Thus, the noise spectra differ 
primarily in the low frequency region where the larger scales of the 
separated turbulence contribute more strongly. 



Figure 4.10: Tail Rotor Harmonic/Broadband Noise Due to Interactions 

with the Fuselage Wake. 
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The engine exhaust flow interaction effects are shown in Figure 
4.11. Both the harmonic and broadband noise are important. The 
effect of turbulent spectral shape (i.e. Dryden vs. Karman) is 
evident. Figure 4.12 shows calculations for different relative 
positions of the tail rotor and the assumed engine exhaust. The 
harmonics are higher for the case where the advancing blade is 
interacting with the maximum velocity increment. 



FREQUENCY (Hz) 

Figure 4.11: Tail Rotor Harmonic/Broadband Noise Due to Interactions 

with the Engine Exhaust Flow. 

The main rotor hub/shaft flow interaction effects are shown in 
Figure 4.13, the wakes of the main rotor hub and shaft are both shown 


to be important. 







Chapter V 


TAIL ROTOR BLADE -VORTEX INTERACTION 

In the preceding chapter, we pointed out that under most 
operating conditions, tail rotors contribute significantly to overall 
helicopter noise in both harmonic and broadband forms. In this 
chapter, we study the tail rotor harmonic noise due to interactions 
between tail rotor blades and the main rotor tip vortices. It is 
found that this noise source accounts for previously not understood 
harmonic component of the noise spectrum. 

Figure 5.1 sketches the geometries of the tail rotor blade-vortex 
interaction and the main rotor blade-vortex interaction. These 
interactions are similar in a sense since, in both cases, the rotor 
blades are interacting with the vortices shed from the main rotor . 
However, they are different because of the relative orientations. The 
tail rotor blade-vortex interactions generally occur with vortices 
with near-normal orientations to the tail rotor plane. 

In the present study, the acoustics of the tail rotor 
blade-vortex interaction is treated by the analysis of Amiet 
( Schlinker and Amiet, 1983, 1986a). The acoustic analysis requires 
information on vortex properties and individual vortex geometry as its 
inputs. In the present study, these inputs are calculated using the 
CAMRAD main rotor free wake geometry analysis (Johnson, 1980; Scully, 
1975) . 
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a) Illustration of Main Rotor and Tail Rotor Blade-Vortex 
Interactions . 



b) Illustration of Main Rotor Blade-Vortex Interaction with Vortex 
Parallel to Leading Edge. 



c) Illustration of Tail Rotor Blade-Vortex Interaction. 

Figure 5.1: Geometry of the Helicopter Tail Rotor Blade-Vortex 

Interactions, from Schlinker and Amiet (1983). 
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The acoustic pressure histories are calculated for selected cases 
for representative helicopters; the results are compared to the tail 
rotor loading and high speed thickness noise calculated by Farassat's 
linear acoustic analysis (Farassat, 1983), and are found to be very 
important. This noise mechanism will be seen to depend strongly upon 
the helicopter operating conditions and the location of the tail rotor 
relative to the main rotor . 

5 i ACOUSTICS OF WEAR-NORMAL BLADE -VORTEX INTERAC TION 

In the present study, the tail rotor blade-vortex interaction is 
modelled as a 2-dimensional blade-vortex interaction problem (an 
airfoil of infinite span chopping through a moving vortex with near 
normal orientation). The analysis was originally developed by Amiet 
(Schl inker and Amiet, 1983; Amiet, 1986a), who assumed that the noise 
is generated by an unsteady loading fluctuation (i.e. dipole) on a 
tail rotor blade when it chops through a near normally incident 
vortex. The unsteady loadings can be calculated using thin- airfoil 
theory, and the far- field noise can then be calculated using Lowson's 
moving dipole theory (section 2.4). 

511 Unsteady Loadings and Far Fi eld Sound 

Since the incident vortices can be viewed as an unsteady upwash 
gust field, we start by treating the general problem of an airfoil 
interacting with an arbitrary gust field. Assume that the flow field 
is linear, which means that the incident unsteady upwash remains fixed 
relative to the mean flow at infinity. After being Fourier 
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decomposed, the z component of the unsteady upwash field in the x-y 
plane can be written in terms of sinusoidal gusts of the form 

w g (x,y) - «( k x .k y ) e ' i(k x (x ‘ Ut)+k y y) (5 L) 

Consider a flat plate of infinite span lying in the x-y plane, 
with leading edge at y - -b (b is half chord of the airfoil). The 
pressure jump across the plate due to the incidence of sinusoidal 
gusts defined by (5.1) can be written as 

Ap(x 1>yi ,k x ,k y ) - 2* Pq U w(k x ,k y ) g(x 1 ,k x ,k y ) e- i <Vl’ k * Ut > (5.2) 

where g( x ^. k x . k y ) is the airfoil gust response function which will be 
given in the next section. 

Using the result of Lowson for a moving dipole (equation 2.29), 
the far field pressure p induced by a point dipole F e^ wt at 
coordinates (x-^,y^,0) can be written as 

p(x,y,z , t) e i(«C+/*(Mx-cO) e-i/iCMx.-Ux.+^yy.)/*) 

4*c a 1 (5.3) 

0 

where 

A* - k x M/yS 2 

a 2 - i-M 2 

a 2 ■ x 2 +/ 3 2 (y 2 +z 2 ) 


The dipole (force) is assumed to be aligned with the positive z 
direction, and the source dimension is assumed to be small compared to 
the source - observer distance. 
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The far field pressure due to the unsteady gust component can 
then be found simply by replacing the force in (5.3) with Ap from 
(5.2), noting that 


- " k x U 


(5.4) 


After integrating over the entire airfoil surface, we have 

o 

in-p U bk x z 

p(k x ,k ) 2 «(k y -^ 2 y/^) w( k x - k y> L(k x ,k y ,M) 

V 

(5.5) 

i(k Ut+/i(Mx-a) ) 

• 6 X 


where 


+b 


L(k x . ky .«) - 4- J 6<x l .^.k y .M) d Xl 


(5.6) 


L is actually the lift per unit span on the airfoil. The exponential 
phase factor in (5.6) accounts for the retarded time effects for the 
sound propagating from source to observer. 

To find the far field pressure - time response, we have to 
integrate (5.5) over all k x and k y gust components. The result is 
given by Amiet (Schlinker and Amiet, 1983) as 


i rbzp U 


2 +co 


p(x,y,z,t) - - 


V J ‘ lk * =(k x'V L < k *'V M) 


CO _ 

0 *® 


(5.7) 


. e i( k x Ut+|i(Mx-a)) dk _ 


where 
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K y ■ »0 2 y /° 

The factor of exp(i^(Mx-cr) ) represents the time At for the sound wave 
to propagate from the airfoil leading edge to the far field observer 
(At - (a -Mx) / /9 2 Cq ) . 

5.1.2 Airfoil Gust Response 

Amiet has shown that the trace velocity (V f ) of the gust for the 
present case is always supersonic (Schlinker and Amiet, 1983). Thus 
the surface pressure gust response function g(^,k x> k y ,M) for a skewed 
gust can be written in terms of that for a parallel gust in 
compressible flow (a 2 -dimensional flow field) using Graham’s 
similarity rule (Graham, 1970). The relation is 

2 

g(/J,k x ,k y ,M) — g(M,k xa) ,0,M oo ) e i(/ik y/ k x^ (5.8) 

where 

k X oo " & 

P 2 - l-M 2 

- H ( l -/ 2 ) 1/2 

* - Mk x //9k y 

^ - M»k x /^ 2 

According to Amiet (Schlinker and Amiet, 1983), the response 
function is divided into two ranges in which two different approximate 
solutions are used. For ^ < 0.4, the response function can be found 
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using a low frequency approximation. The final result for g is given 
by Amiet as 


g(x,k x ,k y ,M) - ±- 



S(k*) e 10** + <f(MJ) 


(5.9) 


where 

• w 2 

f(M) - (l-0)ln(M)+01n(l+0)-ln(2) 


S is the Sears function (Sears, 1940) used in the incompressible 
theory, it can be written as 


S(k) - 


2 

7rk(J 0 (k)-Y 1 (k)-i(J 1 (k)+Y 0 (k))) 


(5.10) 


Yq and Y^ are the Bessel functions of the second kind, subscripts 
denote their orders. With equation (5.6), L can be found to be 


L(k x ,k y ,M) 


\ <> 




(jQ(/jbx/o) -iJ^(M b x/a) ) 


(5.11) 


where J Q and J l are Bessel functions of the first kind, and the 
subscripts denote their orders . 

For the high frequency case where > 0.4, The first two terms 


of g are (Amiet, 1976) 
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M) — e " i(k x (M co( b+x )- 

7r(^k x (l+M oo ) (b+x)) 1//2 

g 2 (x,k k M) - e -i(k x (M co (b+x) -M 2 x 

7r(27rk x b(l+M co ) ) 1//2 

• ( (l+i)E*(2K*M oo (b-x) ) - 1) 


And L can be found to be 


L l< k x- k y' M > 


L 2 (k x ,k yl M) * 


>V2 


Yn E (2e x ) e 

* (k x b(i+M tc )e 1 ) 1 / 2 

(■ 


ie. 


e i0 2 


7r0 1 (27rk x b(l+M a) ) 1 / 2 
2 1 / 2 e- i26 l 


i(l-e' 129 l) + (i 


-* /0 . , M , x 

i n~ E (2pb( ZT + ~ 
(1+(M/M Ux/o)) 1 / 2 M ® r 


where 


0 L - /ib(M m /M - x/a) 

© 2 ■ k x b + /ib(M-x/a) - tt/4 


E is a combination of Fresnel integrals defined by 


E*(x) 


1 x e ’ 

(2,) 1 / 2 0 


dC 


5.1.3 Vortex Velocity Field 


M 2 x-b)/b-7r/4) 

-b)/b-7r/4) 

(5.12) 


r 

-1) E*(4k*M eo ) 
)S> " (5.13) 


Figure 5.2 sketches the geometry of a near-normal blade-vortex 
interaction. The tangential velocity component in the z direction for 
a vortex tilted by 0^ from the z axis whose <j> ** 0 can be written as 
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Figure 5.2: Sketch of 2 -Dimensional Near-Normal Blade-Vortex 

Interaction. 

W ” ~R~ Sin *v (5.14) 

It should be noted that R is the normal distance to the center of the 
vortex, i.e. 

- (xcos^ v - zsin0 v )2 + ^2 (5.15, 

The upwash velocity field due to the vortex can also be expressed in 
terms of the sinusoidal gusts “(k x> K y ). This can be done by taking a 
double Fourier transform of equation (5.14) with respect to x and y: 
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W(k x’ k y> 


sin0. 


( 2 * )' 




II R V S (R) e 1 ( V + V ) dxdy 


(5.16) 


Several vortex models may be used to define the tangential 

velocity field v,(R). The vortex model Amiet used in his analysis is 
given by 


< x + +> r v » (1 ■ e '“ (r/r ° )2 ) 


(5.17) 


where a - 1.25643, r Q v 0 - 172 *, and T is the strength of the vortex. 
Substituting (5.17) into (5.16) and performing the integrals, we can 

get 


w(k x’ k y> 


iky(l+l/2a) tan0 v 


2ir(k^+k^/cos^6 ) 
y x v' 


r 0 v 0 e 


- (k^+k^/cos^^ )r^/4a 
y x v' ' 


(5.18) 


In order to be consistent with the aerodynamic analysis (i.e. 
CAMRAD) , we used a vortex model different from the original Amiet 
analysis. In the present analysis, the tangential velocity v $ is 
defined by the widely used turbulent vortex model 


v. 


2nr 


2 2 
r +r 

c 


(5.19) 


where r c is the vortex core radius (r c is taken to be 0.0025 of the 
main rotor tip radius in the present study). In Figure 5.3, this 
vortex model Is compared to the one used by Amiet. 
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Figure 5.3: Comparison of the Tangential Velocity for Different 

Vortex Models, v^ is normalized by T/4irr c and r is 
normalized by r . 


At large radial distances, the vortex model used in the present 
study decays more slowly than the model used by Amiet. Since a 
different vortex model is used, w, the Fourier decomposed vortex 
velocity field (w considers only the velocity component that is normal 
to the tail rotor disk) , is replaced by 


ik y rr c tan* v 


2 . ,.2 


1 1/2 N 


w(k x ,k ) / 2 2 l/2 - Ki(r c (k; + kVcos%) 7 ) (5.20) 

X Y 4* 2 (k 2 +kVcos Z 9 ) i/Z y 


where K-, is the modified Bessel function of the second kind, and 8 V is 


defined in Figure 5.2. 
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The acoustic pressure - time history for a given blade-vortex 
interaction, calculated by the present vortex model, is compared to 
that obtained using Amiet's original analysis. The comparison is 
shown in Figure 5.4. The results for the two different vortex models 
show close similarities and only minor differences. 



Figure 5.4: Comparison of the Noise Generated from Blade-Vortex 

Interactions with Different Vortex Models. 


In general, <£ v is not zero. Let £ , rj , f represent the coordinate 
system which rotate the original system by <f> v about the z axis. Then 
x, y, z are related to f , r? , ^ by 
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x - | cos - r) sin$ 

v v (5. 21] 

y = £ sin^ v + rj cos<f> v 

z “ f 

Amiet (Schlinker and Amiet, 1983) has shown that in order to 
accommodate various vortex orientations, k x and k y should be replaced 
by k^ and k , where k,- and k are given by 


” k x cos K + k y sin ^v 
k^ - -k x sin<^ v + k y cos0 v 


(5.22) 


It should also be noted that w only accounts for the effect of 
the tangential velocity of the vortices; the axial flow in the main 
rotor tip vortices is neglected in the present analysis, although it 
can be very important as noted in the conclusion. 


52 HELICOPTER MAIN ROTOR TIP VORTEX-TAIL ROTOR INTERACTION 

As noted previously, the acoustic analysis requires inputs from 
an aerodynamic analysis; particularly the information on how the 
vortices are moving on the tail rotor disk, what the properties of the 
vortices are, and what the relative geometry of the vortices and the 
tail rotor blades is. In this section, we will discuss in detail how 
these aerodynamic inputs are obtained. 
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5 - 21 Main Rotor Tip Vortex Free Wake Geometry 

Since the main rotor tip vortex system is generally highly 
distorted under most helicopter operations, classical rigid wake 
analysis can not predict accurate trajectories of the main rotor tip 
vortices (Scully, 1975). The calculation of the free wake geometry of 
main rotor tip vortices is, therefore, very important to the present 
study since the trajectories of vortices directly affect the 
characteristics of the blade-vortex interaction and thus, the noise 
generated. 

In the present study, we use the comprehensive analysis of 
rotorcraft aerodynamics and dynamics code (CAMRAD) by Johnson (1980) 
to obtain information on the properties of the main rotor tip vortices 
and their individual geometries. The free wake geometry analysis of 
CAMRAD is based on a rotor, free wake geometry model of Scully (1975). 

In the present application, the free wake geometry analysis is 
driven by calculated main rotor airloads assuming non-uniform inflow 
at the main rotor disk. It is also assumed that the presence of a 
tail rotor has no effect on the upstream main rotor loading and its 
tip vortex system. Also no fuselage wake or ground effect is 
considered when calculating the free wake geometry. 

5.2.2 Definition of Tail Rotor Blade -Vortex Interactions 

The tail rotor blade-vortex interaction can be calculated 
provided that the following information is given: the normal incident 
velocity of the ingested vortex relative to the tail rotor blade (U) , 
the strength of the ingested vortex element (T) , the skew angle 
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between the ingested vortex element and the tail rotor axis ($ v ) • and 
the ingested vortex orientation (</> v ) . Refer to Figure 5.2 for the 
definitions of U, 6 v , and 4 > v . It should be noticed that these 
parameters are generally not constant as a vortex sweeps through the 
tail rotor disk. 

Since the tail rotor RPM is generally not an integer multiplier 
of the main rotor RPM, the locations of the blade-vortex interactions 
are different for each main rotor revolution. In the present study, 
both #1 blades of the main and the tail rotors are set such that both 
blades will start from - 0° initially. (Figure 5.5 shows the 
definitions of azimuthal angles V> for both main and tail rotors.) The 
information defining a series of blade-vortex interactions can then be 
determined numerically by interpolating the results of the main rotor 
free wake geometry analysis. 

5.3 RESULTS AND DISCUSSIONS 

In this section, we will present several examples calculated for 
UH-1D helicopters to demonstrate the general procedures developed in 
the present study. 

5.3.1 Results on Free Wake Geometry Analysis 

Our simulation begins with the calculation of the main rotor, 
free wake geometry. Results for two speeds are presented; they 
correspond to a UH-1D in level flight at 100 and 60 Knots. The free 
wake geometry results for the 100 Knot case are presented in Figure 
5.6. From the results shown, the interactions between tail rotor 
blades and main rotor tip vortices are evident. 
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5.5. Db f ini t ions of Main Rotor and Tail Rotor Coordinates and 
Azimuthal Angles. 

After the main rotor tip vortex free wake geometry is found, we 
can plot its trajectory/trajectories on the tail rotor disk. Figure 
5.7 shows the possible patterns of the vortex trajectories on the tail 
rotor disk. Figures 5,8 and 5.9 present the actual main rotor tip 
vortex trajectories on the tail rotor disk for our UH-1D examples (100 
and 60 Knot cases). The points shown are interpolated from the free 
wake geometry analysis, and each point represents the location of the 
main rotor tip vortex at an interval of 15° rotation (main rotor 
azimuthal angle) . 
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Figure 5.6: UH-1D Free Wake Geometry Analysis Results, 100 Knot Level 

Flight, Main Rotor V 1 “ 30 . 

Notice that the tip vortices involved in the interactions with 
the tail rotor are shed from main rotor blades at approximately 0 of 
azimuthal position. Also they are relatively "young" (vortex age are 
generally less than 180° for both cases considered), which implies 
that the ingested vortices are not fully rolled-up. (Johnson (1980) 
has suggested that a vortex is not fully rolled-up unless the vortex 
age is larger than 180° or so.) Since a vortex is not fully 
rolled-up, the strength of the ingested vortex should be less than the 
maximum bound circulation on the main rotor blade; we followed the 
assumption made by Scully (1975), and set the strength of the tip 
vortex strength to 0.8 of the maximum bound circulation on the main 


Figure 5.6a: (Continued) Main Rotor ^ - 60°. 


rotor blade span. It should be noted that interactions with vortices 
originating from other azimuthal angles or with some relatively "old” 
vortices are still possible. Such interactions had been observed in 
our calculations for the model scale UH-1, and are certainly possible 
for full scale helicopters operated in appropriate conditions. 

Also the results indicate the lack of multiple interactions, 
i.e., no two vortex traces are found simultaneously on the tail rotor 
disk. This phenomenon has been observed in the model scale experiment 
by Balcerak (1976). But this phenomenon is again possible for full 
scale helicopters, providing the operating conditions and tail rotor 


location are suitable. 
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Figure 5.6b: (Continued) Main Rotor r/> - 90°. 

Our results also show that a certain vortex may be chopped by 
several tail rotor blades. It is quite possible that a certain vortex 
may "breakdown" locally after the first blade-vortex interaction, thus 
a subsequent interaction with the same vortex may be weaker. However, 
the free wake geometry calculations indicate that the main rotor tip 
vortex will generally be drifting through the tail rotor disk. 
Therefore, the subsequent interactions are actually the results of 
tail rotor blades chopping through a fresh piece of vortex. Thus, no 
effects of possible interactions with burst vortices have been 
included in the present study. 


Vortex 


Vortex 1 


Vortex 2 



(a) Single Vortex Interaction 



( c ) Two Vortex Interactions - 
Two Trajectories 


(b) Two Vortex Interactions - 
Single Trajectory 



(d) Complex Interaction - 

Multiple Vortex/Trajectories 


Figure 5.7: Possible Main Rotor Tip Vortex Trajectory Patterns on the 

Tail Rotor Disk. 


5.3.2 Acoustic Signals from the Blade Vortex Interactions 

Figure 5.10 shows the results for the tail rotor blade-main rotor 
tip vortex interaction noise signals for the 100 Knot casej horizontal 
tick marks in the figure are 0.1 seconds apart. Figure 5.11 shows the 
same interactions as in Figure 5.10 with a higher time resolution 
(tick marks are 0.01 seconds apart). Figures 5.12 and 5.13 show 
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Figure 5.8: Main Rotor Tip Vortex Trajectory on Tail Rotor Disk, 

UH-1D Helicopter, 100 Knot Level Flight. 



Figure 5.9: Main Rotor Tip Vortex Trajectory on Tail Rotor Disk, 

UH-1D Helicopter, 60 Knot Level Flight. 
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similar results for the 60 knot level flight case. One important 
aspect of the tail rotor blade-vortex interaction we see from these 
results is that the strength of the acoustic signal does not decrease 
when the helicopter is flying at slower speed. Moreover, the 
interaction seems to be more frequent in such case. It should be 
noted that the above calculations assume a bottom- forward rotating 
tail rotor. 



Figure 5.10: Sound Pressure History for Tail Rotor Blade Vortex 

Interactions, UH-1D, 100 Knots. 

5-3.3 Comparison with Other Harmonic Noise Mechanisms 

To determine the relative importance of the tail rotor blade-main 
rotor tip vortex interaction, results are compared to the tail rotor 
loading and to high speed thickness noises. The loading and thickness 
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Figure 5.13: Sound Pressure History for Tail Rotor Blade Vortex 

Interactions, UH-1D, 60 Knots. 

noises shown are obtained with a discrete frequency noise prediction 
program WOPWOP using Farassat's linear acoustic analysis (Farassat et 
al., 1980, 1983). The airloads on the tail rotor blades are obtained 
with approximate aerodynamic analysis. Figures 5.14 and 5.15 show the 
comparisons of loading and thickness noise with some representative 
blade-vortex interaction signals for the 100 knot case. Clearly the 
blade-vortex interaction signals are very important. These figures 
also show the detailed shape of a typical tail rotor blade-main rotor 
tip vortex interaction. 
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5.3.4 Effects of Changing the Tail Rotor Location 

As discussed previously, the vortex trajectory on the tail rotor 
disk is very important to the tail rotor blade-vortex interaction 
noise. The tail rotor location relative to the main rotor, and the 
helicopter operating conditions are two primary factors that determine 
the vortex trajectories on the tail rotor disk. For the cases shown 
so far, the interactions usually occurred at the retreating side of 
the tail rotor; this effectively reduces the incident velocity U, and 
therefore the strength of the interaction noise. To study the effect 
of tail rotor location on the blade-vortex interaction noise, we 
artificially lowered the UH-1D tail rotor by 0.5 m. This will cause 
the blade -vortex interactions to occur on the advancing side of the 
tail rotor, thus enhancing the strength of the interactions. 



Figure 5.16: Main Rotor Tip Vortex Trajectory on Tail Rotor Disk, 

UH-1D Helicopter, 100 Knot Level Flight, Tail Rotor 
Lowered by 0 . 5 meters. 
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For the 100 knot level flight case, the main rotor tip vortex 
trajectory on the tail rotor disk is now shown in Figure 5.16. Notice 
that the path is higher than that shown in Figure 5.8 due to the 
lowered tail rotor. The acoustic pressure - time history of the 
resulting blade-vortex interaction is shown in Figure 5.17. Notice 
that there are considerable differences between the results shown in 
Figures 5.10 and 5.17. Since the vortex is passing through the 
advancing side of the tail rotor, the relative velocity between the 
tail rotor blade and the vortex element is generally higher , so the 
interaction signals have higher magnitude. Also, the interactions are 
more frequent than in the previous cases. Unquestionably, with this 
configuration (with lowered tail rotor) , tail rotor noise would be 
higher than that from a standard tail rotor. 

5.3.5 Comparison with Experiments 

As mentioned previously, very few experiments have ever addressed 
the tail rotor noise problem. The most important experiment we have 
examined is that of Balcerak (1976). In his parametric studies of 
tail rotor noise mechanisms, he found that the noise increased as the 
main rotor wake passed the tail rotor disk. In the present study, we 
examine several of his cases. Presented here is a special case for 
which the tail rotor RPM is an integer multiplier of the main rotor 
RPM. In this case, the interaction will be periodic. The 
corresponding advance ratio is 0.2, and the model is an l/16th scaled 


UH-1 series helicopter. 
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Figure 5.17: Sound Pressure History for Tail Rotor Blade Vortex 

Interaction, UH-1D, 100 Knots, Tail Rotor Lowered by 
0.5 meters. 


Figure 5.18 shows the free wake geometry analysis results; unlike 
the full scale cases, multiple interactions (with vortices originating 
from approximately 0 and 180°) do exist in this example as mentioned 
by Balcerak in his study. The 180° vortices are also quite "old" 
(their wake age are about 540°) . Figure 5.19 shows the tip vortex 
trajectories on the tail rotor disk. The acoustic signals due to the 
blade -vortex interactions are shown in Figure 5.20. Figure 5.21 shows 
the acoustic pressure spectrum obtained from the above results along 
with the experimental results of Balcerak. (Notice that the 
calculated results shown do not include the loading and high speed 
thickness noise.) 




Figure 5.18: Model UH-1 Free Wake Geometry Analysis Results, Level 

Flight (n - 0.2), Main Rotor V» “ 30 . 




Figure 5.18a: (Continued) Main Rotor 0 - 60 . 
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Figure 5.20: Sound Pressure History for Tail Rotor Blade Vortex 

Interaction, Model UH-1, /i — 0.2. 



Figure 5.21: Comparison of Tail Rotor Blade Vortex Interaction Noise 

Spectrum and the Experimental Result of Balcerak (1976) 
Model UH-1, u = 0.2. 



Chapter VI 


CONCLUSIONS 

The understanding of and the ability to predict broadband rotor 
noise are approaching a satisfactory state in many respects. The 
present study shows that the important broadband noise mechanisms are 
now understood well enough to be able to make predictions to within 
about 5 dB of experimental data. This understanding should enable 
designers to minimize broadband noise in the cases where it is a 
controlling factor in a design. 

The calculations and comparisons made in the present study 
indicate that inflow turbulence induced lift fluctuations are the most 
important broadband noise sources at low frequencies. This radiation 
can be predicted down to the lowest blade passing frequency, including 
the smooth peaked spectral structure, by the method of Homicz and 
George (1974). For the higher frequencies, which are of more 
practical interest, the methods of George and Kim (1976) and of Amiet 
(1976) are just as satisfactory and are much easier to compute the 
method of Homicz and George. When the same inflow turbulence spectrum 
is used, both methods agree well with each other and with measurements 
over a full range of parameters, except at angles within about 10° to 
rotor plane where other in-plane mechanisms are dominant. 

The Karman spectrum, which has been implemented in both George 
and Kim's and Amiet' s methods, is shown to be suitable for use in 
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predicting the inflow turbulence noise from full-size rotors. 

However, the Dryden spectrum, which is only available in George and 
Kim's method, is shown to be more suitable to predict the inflow 
turbulence noise for indoor model rotor where small scale turbulence 
is involved. The accuracy of the prediction depends strongly upon the 
description of the inflow turbulence: the integral length scale, A, 

the turbulent intensity, w, and the turbulent spectral shape. Further 
research on turbulence properties is needed in the future. 

Boundary layer trailing edge noise is now well understood. The 
analyses of Kim and George (1980) and of Schlinker and Amiet (1980) 
and the correlation of Fink (1978) all appeared to give results which 
agree reasonably well with experiments. This source often is the most 
important noise source at high frequencies on large rotors when inflow 
turbulence is weak. It increases significantly with angle of attack 
due to the increase of boundary layer thickness. 

Tip vortex formation noise is satisfactorily predicted, although 
it does not dominate in any of the experiments to definitively 
establish the precise accuracy of the model of George et al . (1980, 

1983b, 1984a). More experimental data is needed on flows and 
fluctuating pressure on different shapes of rotors and wing tips. Tip 
vortex noise is most important for rotors with square tips and wide 
chords at high angles of attack. This noise can be reduced 
significantly by detailed tip shape changes, but this is presently 
unexplored . 

Turbulent vortex shedding noise from blunt trailing edges is now 
reasonably predicted. Generally speaking, the noise spectra due to 
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turbulent vortex shedding from blunt trailing edges show a peak 
occurring at certain frequency range. The peak frequency depends upon 
the trailing edge thickness; a small trailing edge thickness will 
generate a high frequency peak and a thick trailing edge will result 
in a peak of lower frequencies. The level of the spectrum peak also 
depends on the thickness of trailing edge; peak level increases 
roughly to the third power of trailing edge thickness, thus a slightly 
blunted rotor trailing edge can contribute significantly to the 
overall noise. More measurements of surface pressure fluctuations 
near blunt trailing edges are needed to improve the accuracy of the 
prediction. 

helicopter ta.il rotor broadband noise, a number of 
^PP ima te analyses of the interaction noise sources bave been 
presented. The interaction noises are generally much greater than the 
isolated tail rotor noise. Inflow turbulence noise is the most 
important tail rotor broadband noise mechanism: the turbulent wakes of 
the main rotor and the fuselage are the most important sources of the 
tail rotor turbulent inflow, while the turbulent main rotor tip vortex 
is of less significance. Further studies are needed to determine the 
turbulent properties in the tail rotor inflow. 

The present study has shown that fuselage, hub, and hub shaft 
wakes and engine exhaust may be important to harmonic noise if the 
tail rotor is fully immersed in the respective disturbed flows. In 
order to minimize noise, it is desirable to place the tail rotor where 
it will not be affected by any of the upstream disturbances. While 
the engine, fuselage, etc. flows may be avoided or reduced, the main 
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rotor wake is so large that it tends to interact the tail rotor over a 
wide range of flight conditions. The flow disturbances in the main 
rotor wake (including both the overall downwash field of the main 
rotor mean wake and the main rotor tip vortex wake) are found to be 
the strongest tail rotor harmonic noise sources. 

The interaction between the tail rotor blade and main rotor tip 
vortices is very important to harmonic noise. The noise generated 
depends strongly upon the helicopter operating conditions and on the 
tail rotor location. Major parameters governing this blade-vortex 
noise generation are the ingested vortex strength, the ingested vortex 
orientation relative to the blade, and the relative velocity of the 
ingested vortex to the tail rotor blade. The present study has shown 
that this noise mechanism is at least of the same order of magnitude 
as some of the strongest tail rotor noise sources such as high speed 
thickness noise. Future studies should consider a vortex chopped by 
an airfoil of finite span and the possibly major effect of the axial 
flow in the main rotor tip vortex. 
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